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Abstract 



We try to review the main current ideas and points of view on the 
running coupling constant in QCD. We begin by recalling briefly the 
classic analysis based on the Renormalization Group Equations with 
some emphasis on the exact solutions for a given number of loops, in 
comparison with the usual approximate expressions. We give partic- 
ular attention to the problem of eliminating the unphysical Landau 
singularities, and of defining a coupling that remains significant at 
the infrared scales. We consider various proposal of couplings directly 
related to the quark-antiquark potential or to other physical quanti- 
ties (effective charges) and discuss optimization in the choice of the 
scale parameter and of the renormalization scheme. Our main focus 
is, however, on dispersive methods, their application, their relation 
with non-perturbative effects. We try also to summarize the main 
results obtained by Lattice simulations in particular in various MOM 
schemes. We conclude briefly recalling the traditional comparison with 
the experimental data. 
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1 Introduction 



The renormalization procedure in field theory consists in a redefinition of the 
unrenormalized constants that appear in the Lagrangian, in such a way that 
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the observable quantities can be kept finite when the ultraviolet cut off Auv 
is removed. In the framework of any given renormalization scheme (RS), 
dimensional reasons require the introduction of a new quantity fi with the 
dimension of a mass. Note that observables must be independent of the par- 
ticular RS and of fi by definition. Intermediate quantities, like renormalized 
coupling constants, masses and field normalization factors depend on RS and 
on fi by construction and a change in RS amounts to a redefinition of such 
quantities. Obviously, approximate expressions of the observables depend in 
general on u, and RS. 

In QED the masses of the charged particles have a direct physical meaning, 
so there exists a natural scale referring to which fi can be fixed. In QCD, 
due to confinement, such a natural scale does not exist and RS and fi have 
to be chosen with other criteria, those of simplicity and of convergence. 
In QCD we have a single coupling constant g s , or the usually more convenient 
a s = j^, and various quark masses mj with f — u, d, . . . , t. We refer to 
their dependence on fi in the framework of a given RS (a s (/i 2 ), mj(/i 2 ), . . .) 
as to the running coupling constant, to the running masses and so on. MS 
is the simplest and most commonly used scheme but, being essentially per- 
turbative, alternative definitions of a s (fi 2 ) may be more appropriate in many 
cases. 

Even in QED perturbative expansion is believed to be only asymptotic. As 
known this means that the approximation to the considered quantity im- 
proves as the number of terms included increases, until a maximum number 
Af* is reached. After this the single terms become progressively larger and 
the series loses any meaning. In QED a ~ 1/137 and Af* ~ 1/a ~ 137, 
so in practice no problem arises from this lack of convergence. However in 
QCD a s is nearly two orders of magnitude larger and AT* is of order 1. An 
appropriate choice of fi and possibly of RS becomes therefore essential. 
The fi dependence of the renormalized quantities is controlled by the renor- 
malization group (RG) equations. Let us concentrate on coupling constant 
a s (fi 2 ) which is related to the unrenormalized constant, a" r , by the usual 
equation 



V Jv uv/ 

where A^y is the cut-off parameter. The corresponding differential equation 
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dfi 2 
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can be obtained differentiating eq.flf . /3(a s ) remains obviously finite as 
Auv ~~ > oo and in perturbation theory takes the form 

(3{a s ) = -a 2 {(3 + fta B + (3 2 a 2 + ...). (3) 

As is known, the various terms in /3 , j3\, /3 2 , ... correspond to one loop, 
two loops, three loops . . . contributions; (3q and (3\ are universal in the mass 
independent schemes. 

Note that, for a general RS, Z ai /3(a s ), /3o, Pi, ■ ■ ■ depend also on quark 
masses through the variables rrij/fj?, not explicitly indicated. However, ac- 
cording to the decoupling theorem, all quarks with masses much larger than 
the energy scale of interest (in particular rrif ^> fi) can be ignored. On the 
contrary, if rrif <C fi, we can often neglect rrif. Then, the discussion can be 
greatly simplified if for every fi we divide the quarks in active quarks with 
rrif = and inactive ones, which we simply ignore. Within this framework 
Po, /3i, ... depend on /i only through the number of active quarks rif, which 
changes by ±1 any time ii crosses a quark threshold rrif . Furthermore the 
first two coefficients, (3q and are RS independent, while all the others 
depend on the scheme. In the one loop approximation (i.e. keeping only the 
first term in (jSJ)) eq.© gives 

1 + /? as(/io) m (/^7^o) ^ V fit J 

which explicitly expresses a s at the \x scale as a function of the same quantity 
at the fj, scale. Eq. ^ clearly shows that a change in the value of \x consists 
in a reorganization of the perturbative expansion of any observable or, what 
is the same thing, in a resummation of various contributions. Setting 



A 2 = /igGXP 



f3 a s (n 2 c 



(5) 



ce s (fi 2 ) can be written in terms of the overall scale A, without any reference 
to a specific /i 

as(/i2) = /3 ln(ivA 2 )" (6) 
As concerns the best choice of fi 2 in a specific calculation, let us consider the 
perturbative expansion of an amplitude or observable G(q, x) of canonical 
dimension 0. We assume G written in terms of an overall momentum q and 
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a set of dimensionless variables (angles, Bjorken variables and so on) which 
we collectively denote by x. According to the process we are considering q 
may be space-like, q 2 < (a momentum transfer), or time-like, q 2 > (an 
energy). We shall often set q 2 = =fQ 2 (Q 2 being as a rule positive) or also 
q 2 = s in the time-like case. 

Under the above mentioned assumption that the active quark masses can be 
neglected, we can write 

As we mentioned, G must be independent of /i 2 , so 

Then by replacing (J7J) and (J3J) in (JHJ), we obtain a system of differential 
equation in the coefficients of (J7J), whose solution gives (see e.g. sec. 3.3 for 
details) 

#o = 9 (x) 

92 = 9 2 (x) 

93 = g 3 (x) 

Let us assume that g (x), ~g~i(x), ... decrease sufficiently fast in order ((7|) to 
become significant when Q 2 ~ u 2 . This would clearly no longer be true for 
very different values of Q 2 . We must obviously modify fi 2 choosing it always 
on the scale of interest. In particular it may be convenient to chose exactly 
fi 2 = Q 2 and then (JTJ) takes the form 

G{q,x)=g Q {x)+g 1 [x) + g 2 (x) +..., (10) 

showing that a s (Q 2 ) is actually the most convenient expansion parameter for 
calculating the quantity G(q,x). 

Now let us go back to eq. (JU). We have (3q = ^-(11 — 2n//3) and, since 
the number of presently known flavors is 6 , this remains positive at least 



9i = 9i(?) 

- n(3 g~i(x) hr % 

- (ifPig^x) + 2nf3 g 2 (x)) In % + ^PoU*) ^ % 

ur n, A 



(9) 
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in the entire range so far accessible. Then, for fi — > oo we have a s (fi 2 ) — > 
(asymptotic freedom) and as we can see (see sec. 2) this remains true even 
if a larger number of terms is taken into account in (J3J). 
Conventionally the mass of the Z boson, M z ~ 91.2 GeV, is used as 
Ho reference value. The world average of all determination is currently 
a s (M|) = 0.1189 ± 0.0010 EEE|. 

With diminishing Q 2 we may expect expansion ((7j) or (jlOj) to remain mean- 
ingful up to a few GeV if q is space-like. The situation requires more attention 
for time-like q, q 2 = s > 0, as discussed in sec. 2.5. Then, in fact, in (JHJ) 
(3o In % — > Po In ~ S T° = A) (In -% — Z7r) and, even if we take fx 2 = s to control 

fl fl fl 

the logarithms, now the coefficients g 2 , g 3 ... are modified for terms propor- 
tional to powers of fto^, which become rather large as the order increases (see 
sec. 2.5). Then the "space-like" coupling a s (s) remains a good expansion 
parameter for large s, but for some intermediate s may be more appropriate 
to define a new "time-like" coupling a(s) in which the terms in /3ovr are reab- 
sorbed, so that the coefficients remain small. This can be easily done starting 
from eq. (jlOj) rather than (J7J) and considering the analytic continuation of the 
coupling rather than the coefficients. Naturally using the time-like coupling 
rather than the space-like one can be interpreted as a change of RS and, as 
such, the two expressions become identical for large s. 

Note furthermore that a s (/i 2 ) , as given by (0), has a pole in [i 2 = A 2 . In 
the MS scheme, the expression remains singular even if more loops are taken 
into account in Q, and the singularity structure changes. Such singularities 
are obviously related to lack of convergence of the series occurring in (j3J) but 
they are clearly non physical. In fact in a generic scheme and in a generic 
gauge a s (/i 2 ) is a kind of intermediate quantity only indirectly related to 
actual observables and it is often chosen with criteria of formal simplicity. 
It is just the case for the mentioned MS scheme and for the Landau gauge 
to which as a rule we shall refer. On the contrary we know that on general 
grounds G(q, x) is analytic in the entire q 2 complex plane apart from a cut on 
the real positive axis from some threshold to +oo and in its exact expression 
the above non physical singularities must cancel. Correspondingly, eq. fTUj) 
suggests that it is possible to define schemes in which a s (Q 2 ) can be extended 
again to the entire complex Q 2 plane apart from a cut on the real negative 
axis. To solve the problem of the elimination of the above spurious singu- 
larities various strategies have been proposed. We remind here for instance 
the Brodsky-Lepage-Mackenzie (BLM) criterion 3\, which fixes the scale \x 
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such that the next-to-leading order (NLO) coefficient of perturbative series 
is independent of n/ , and also the so-called "Fastest Apparent Convergence" 
(FAC) technique [I], that amounts to setting the scale \x so that the NLO 
coefficient and all the higher ones are zero. The latter is related to the ef- 
fective charges jl], which consist in defining new coupling constants 
«eff(Q 2 ) in more strict connection with some physical observable (see sec. 
3.2). Typically let us assume that the quantity G in (jTUJl is independent of 
x (or x has been fixed on some special value) and set exactly 



Effective charges are obviously well defined everywhere in principle and easily 
extracted from the data. They depend on the particular observables chosen 
and, accordingly, can be either of the time-like or space-like type. However, 
they can be related to each other by referring back to the MS scheme. In 
fact, comparing (fTTj) and (JTUJ), we have 



Another device is the Optimized Perturbation Theory (OPT) [Z2]-[ZH], which 
consists in improving the convergence of the perturbative expansion by choos- 
ing the value of u 2 and of some of the RS parameters (like (3 2 , P3 themselves) 
at every order according to a criterion of minimum sensitivity. This has been 
discussed as a matter of example in some details in sec. 3.3, with reference 
to the theoretical predictions it yields on the IR behavior of the coupling. 
Finally another method , summarized in sec. 4, consists in requiring ab initio 
the desired analyticity properties for a s (/x 2 ) , by rewriting this as a disper- 
sion relation and applying the perturbative theory to the spectral function 
Im [a B (— m 2 — iO)} rather than to a s (/i 2 ) itself (ZHJ-IHII- In this way Landau 
singularities are suppressed at their very roots. In the one loop approxima- 
tion from eq.(@} we find 



G(Q 2 )=g + g- 



a e s(Q 2 ) 



(11) 



TT 




(12) 



Im [a s (—m 2 — iO)] 



7r/?oO!s(^o) 



(13) 



[1 + A)Q! s (^o) ln ( m2 / '^o)] 2 + (^Otts^o)) 2 



which is finite everywhere on the positive m 2 axis and can be considered in 
some way as a reorganization of the series occurring in the last term in eq.(0} 
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after evaluation of the imaginary part. 

It must be stressed that all these techniques lead to a finite a s (0) and are 
essentially consistent. Furthermore all of them use as input the coefficients 
of the power expansion of a s (n 2 ) itself or of specific observables in terms of 
a s (nl). In this sense we may call them perturbative 1 . 

However, we believe that there are intrinsically non perturbative effects in 
the theory. These are in part related to the Gribov problem of the multi- 
ple solutions of the gauge fixing equation and more seriously with the 
occurence of the singularities in a s = which can not be obtained by any 
manipulation of perturbative expansions. Confinement and string tension, 
condensates, topologically non trivial classical solutions are of this type. In 
this context, unfortunately, we can not rely on any rigorous methods in the 
continuum. We have achieved some intuitions and very partial results, but 
we are far from a comprehensive understanding of the all the matter. Within 
this framework, what we can do is to parametrize such non perturbative ef- 
fects in terms of some universal quantities that occur in different types of 
phenomena. In particular, according to perturbation theory, in high energy 
processes observables should vanish logarithmically with the energy (apart 
from mass corrections). Non perturbative effects are expected to appear in 
terms of powers of the 4^ and ^ InQ 2 form (see sec. 4.5). Such terms should 
in principle be detected at somewhat intermediate energies even if the ex- 
perimental situation is not sufficiently clear. In this context the coupling is 
written as the sum of a perturbative af T (Q 2 ) and a non perturbative part 
af F (Q 2 ). The latter is supposed to vanish fast enough out of the infrared 
region and the coefficients A u , A' u can be expressed in terms of moments of 
af F (Q 2 ) or of its time-like counterpart af p (m 2 ) [Hi]. 

The only really first principle method that we have to handle non perturba- 
tive effects is to approximate the continuum with a lattice and then rely on 
numerical simulations. As applied to the running coupling this method de- 
livers results in substantial agreement with those obtained in the continuum. 
The significance of the latter result is however uncertain, due to finite size 
limitations. 

In this review the emphasis is put on the IR behavior of the coupling con- 

1 However, at the one loop level, the analytic coupling can be written explicitly in terms 
of the ordinary coupling in the form (see later) a an = a s + (l — e f, ° as ^ that has an 
essential singularity in a s = and in this sense can be said non perturbative. On the 
contrary the relation between a an (/^ 2 ) and a an (/j,Q) is analytic. 
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stant and in particular on the analytic methods. We have tried to give a 
unified presentation of the matter. To this end we often took the liberty of 
modifying notations used in the original papers. Due to the extent of the 
literature on the subject, and to space limitation, we had to make appro- 
priate choices. In this sense we do not claim our paper to be exhaustive or 
complete. The choice we have made is obviously due to our personal taste or 
to the specific point of view we have adopted. In particular in some cases we 
have discussed with a greater detail some specific papers as an example of 
a more general methodology. In such a context we apologize in advance for 
giving less attention to some papers than to others, even very significant. As 
a rule, we have kept to a general theoretical level and referred to application 
only occasionally. Apart from a brief summary of significant results in the 
concluding remarks, no attempt was made to discuss the actual methods by 
which the coupling is extracted from experiments at the various energies. For 
this we refer interested people to specific reviews on the subject such as jH] 
and [IE]. 

The plan of the paper is as follows. In sec. 2 we have summarized the main 
results of the usual perturbation theory as reference. For details and bib- 
liography on this subject we refer to standard books and reviews (see e.g. 
P3 12 113 Ell)- We have tried to emphasize the importance of exact solutions 
to renormalization group equations with a given number of loops in compari- 
son with the usual iterative solutions appropriate for high energies. This was 
due to our interest in the IR region 2 . Some attention has also been given to 
the need to define a specific time-like coupling. In sec. 3 we have considered 
some purely phenomenological modifications of the coupling constant expres- 
sion, suggested essentially by the quark- ant iquark potential theory. We have 
also discussed some significant examples of physical couplings, a classic ap- 
plication of the optimized perturbation theory. Sec. 4 is entirely devoted 
to the dispersive approach, analytic coupling constant and relation between 
space-like and time-like expressions, the parametrization of non-perturbative 
effects and the so called Analytic Perturbation Theory. In sec. 5 we have 
tried to summarize the present state of the studies on the running coupling 
on the lattice. In sec. 6 we outline our conclusive remarks and also discuss 
briefly some information on the coupling under 1 GeV that can be obtained 



2 In this connection it should be recalled that the expansion in numbers of loops is an 
expansion in the Planck constant h and is conceptually different, even if related, to a power 
coupling expansion. 
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from relativistic studies of the light quark- ant iquark spectrum and from a 
parametrization of other non-perturbative contributions. 



2 Perturbation Theory 



2.1 QCD p function 

As we mentioned, eq. (J2J) is derived applying the operator fi 2 (d/dfi 2 ) to 
eq. (0). Thus one needs to compute the renormalization factor for the cou- 
pling Z a and this can be accomplished in several ways. One can start from 
the quark-gluon vertex Zq qg together with renormalization factors of quark 



and gluon propagators Z q and Z„ to obtain Z 



y2 7-27-I 



but other 



choices, as ghost-ghost-gluon vertex or trilinear and quartic gluon interac- 
tions, equally work. If dimensional regularization is used the limit Auv ~^ 00 
can be anticipate as long as the space dimension D = 4 — e is kept dif- 
ferent from 4. Taking into account the physical dimension of the fields in 
.D-dimensional space the coupling constant acquires a not vanishing physical 
dimension, thus it is usual the parametrization a™ = fi £ a s Z a , where in the 
MS renormalization scheme 



n=l 



(14) 



and then 



0(a s 



2 s da s a 



(15) 



From this equation and the explicit perturbative form of Z a we have [12 



1 

47T 



Po = — 



01 



02 



03 



(4tt) 2 
1 

1 

(4^)4 



11-3% 



102 -— n, 
3 ; 



2857 5033 325 ' 
s 54 f 



2 1 
149753 

— 4 



3564C; 



1078361 6508 



162 



27 



6 ) 



Hi 



( 50065 6472 \ 2 1093 , 



V 162 



729 



(16) 
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where ( u is the Riemann zeta-function, £3 ~ 1.202057. The coefficients (3j 
generally depend on the RS employed, but the first two are universal among 
the massless schemes. Moreover, in the MS-scheme the /3-function is gauge- 
independent at any order [7j, and in an arbitrary mass- dependent scheme 
this feature is preserved only at the first order. 

As well known, the universal one-loop coefficient jS] has a positive sign pro- 
vided there is a small enough number of quark fields (rif < 33/2), thus the 
theory is asymptotically free, that is the /3-function has a stable UV fixed 
point as its argument approaches zero; indeed this coefficient is the sum of 
two contributions, the relevant one with respect to asymptotic freedom prop- 
erty being the first, which arises from pure gauge field effects i.e. from the 
nonlinear Yang-Mills interaction terms. The two-loop coefficient has been 
computed for the first time in |U] and is positive up to rif = 8 . 
Higher order approximations are scheme- sensitive, and it is common practice 
to perform computation within MS or MS procedures which have the same 
the /3-function. The first calculation of three-loop coefficient is due to [TU] . 
where the ghost-ghost-gluon combination in the Feynman gauge was used. 
In a more recent work ^Tj the quark-gluon vertex was instead employed, pro- 
viding an independent check in an arbitrary covariant gauge of the previous 
result and its gauge-independence. Finally, the original four-loop calculation 
|12'] has been performed using the ghost-ghost-gluon vertex in a arbitrary 
covariant gauge, and for a generic semi-simple compact Lie symmetry group. 
The result turns out to be gauge-independent as expected within MS pro- 
cedure, and involves higher order group invariants such as quartic Casimir 
operators; specialized to the standard SU(3) symmetry, the four-loop coeffi- 
cient is a positive number for every positive rif (see also [T3j ). Finally note 
that all four coefficients are positive up to rif = 6 except for /3 2 which is 
negative for 6 < rif < 40 . 

2.2 Running coupling 

Evolution of QCD running coupling can be gained integrating the differential 
equation (0), that can be rewritten as 



The exact one-loop solution or (JHJ) is obtained by straightforward integra- 
tion retaining only the first term in (jH)). We report for reference the second 




(17) 



12 



form 



ct s (/i 2 ) 



1 



(18) 



Aln^/A 2 )- 



As yet noted the dimensional scale A keeps track of the initial parametriza- 
tion (hq, a s (/ig)), and it is scale-invariant; its value is not predicted by the 
theory but must be extracted from a measurement of a s at a given reference 
scale. Emergence of a scale parameter, sometimes referred to as dimensional 
transmutation, breaks naive scale invariance of the massless theory, and it 
is commonly believed to be associated with the typical hadron size i.e. to 
the energy range where confinement effects set in. Roughly speaking, A is 
the scale at which the (one-loop) coupling diverges (Landau ghost), and per- 
turbation theory becomes meaningless. Further, it is scheme- dependent and 
receives further corrections at each loop level, but for simplicity we use the 
same notation throughout. 

In the next loop level the integration of (|17|) leads to a transcendental equa- 
tion, that is, starting from the two-loop approximation to the /^-function in 
(j3J), straightforward integration in (|T7j) yields 



where B\ = /?i//?q and the constant term from the lower end points can be 
again reabsorbed into the A-parametrization, with the commonly adopted 
prescription (see e.g. [T^ fT5|) 



which fixes a specific choice for A . Thus we get the two-loop implicit solution 



from which the two-loop scaling constant is immediately read with fi = fi^. 
To achieve an explicit expression for the running coupling at this level one 
should resort to the many- valued Lambert function W(() defined by 





In — = C-EilnA 
Mo 



(20) 




(21) 



W(C)exp[W(C)]=C, 



(22) 
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which has an infinite number of branches Wk(C) k — 0, ±1, ±2 . . . such that 
Wn(0 = W^-n(C*) (f° r more details we refer to [IE])- The exact solution to 
eq. (}2*Tj) . being B x positive, for < rif < 8 reads 3 (see e.g. fl"7] ) 

1 1 1 / 1 \ 1//Bl 

/3 0j Bi 1 + tV-i(C) eB 1 \zJ 

where 2 = fi 2 /A 2 , and is the "physical" branch of the Lambert 

function, i.e. it defines a regular real values function for ( G (— e _1 , 0) which 
fulfills the asymptotic freedom constraint. Indeed, W_i(C) as a function of 
complex variable has a branch cut along the negative real axis (actually a 
superimposition of two cuts starting from — oo to — e _1 and to respectively), 
and defining it on the cut coming from the upper complex plane it assumes 
real values in the interval (— e _1 , 0) , with the asymptotics 

W_i(-e) =m£ + 0(m|m£|) (24) 

W-x + ej = -1 - V2ee + 0(e) (25) 

as e — > + . Outside this region of the real axis it takes on complex values, 
indeed it is not a real analytic function. Though not easy for practical aims, 
eq. (J23)) yields the most accurate expression for investigating the IR behavior 
of the running coupling, since it has not been derived by means of deep 
perturbative approximations (aside from the truncation of the two-loop (3 
function). Actually, a frequently used two- loop approximate solution, known 
as the iterative solution ([10]), is obtained starting from eq.(fT9|) together with 
a single iteration of the one- loop formula (|18|) . that is 

«f (*) = ^ 1 r , (26) 

:t V ' lnz + £iln(l + B^\nz) V ; 

where z = /i 2 /A 2 , and A now defined by 

ln^ = C-Siln^ (27) 



o 



3 Note that if9<n/<16 the principal branch Wo is involved, but here and throughout 
our discussion is focused on the physical values n/ < 6 . 
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in (JE2), related to the standard one ({20} by the factor 



ln(A/A) = -5 1 In B 1 . 



(28) 



However, the commonly used two-loop solution is an asymptotic formula 
which strictly relies on the smallness of a s for fairly large fi 2 , since it amounts 
on solving eq. (|19|) (with the choice (|2U|) ) where the last term in the r.h.s. has 
been neglected. Again after one iteration of the one-loop formula, the result 
is then re-expanded in powers of 1/L, where L = \nz and z as before 



a?Hz) 



Po \nz 



Pi \nz 



(29) 



Eq. ()29|) is known as the standard two-loop running coupling and works only 
in the deep UV regime, i.e. for L ^> 1 (see fig. 1(a)). 

Under the same assumptions one can easily derive the three and four-loop 
approximate formulas; at any loop level asymptotic RG solutions are ob- 
tained as a rule through a recursive recipe involving the previous order result 
as an input in the transcendental equation arising from integration of 1)17]) . 
Starting with the approximate implicit solution at four-loop level (obtained 
by expanding the integrand on the r.h.s. of eq. (j!7)l ) 



In — — 7t 

Mo 



G + - h^ln ot s + - 



Pi , Pf-2p p 1 p 2 + p 2 p 3 



p a s p 2 



-« s 



a 2 + 0(a 3 s 



(30) 

the analogous trick which led to eq.(j2Hl) yields the UV asymptotic four- loop 
running coupling in the standard form of an expansion in inverse powers of 
the logarithm L for L ^> 1 (see e.g. |27| ) 



(4) 



P L 



Px InL 

Pl^r 1 p 2 l 2 



+ 



d(ln 2 L-hiL 



PIV 



-ln 3 L + -ln 2 L + 21nL 
2 



0o 
1 

' 2 



M . T , Ps 

3^^1nL + — 



01 



+ 



(31) 



which turns out to be nearly indistinguishable from the three-loop curve 
(see fig. 1(b)). Here the one- loop solution eq.tflHjl has been emphasized, and 
being the leading UV term in (|31)L it defines the asymptotic behavior of 
the perturbative running coupling, i.e. it clearly exhibits the asymptotic 
freedom property. On the other hand, the two and three-loop asymptotic 
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Figure 1: (a) We show the fractional difference between ea. (|23|) and ecis. lf2"§l) . lf2f)|) (solid 
and dashed line respectively) , with A = 350 MeV and rif = 4 . (b) From 6 : it is dis- 
played the fractional difference between the 4-loop and the f-, 2-, 3-loop approximations 
to eq.Jjnj, with n/ = 5 and normalizing conditions for all curves at a s (M§) = 0.119 . 

expressions are easily read from eq. (|31j) by keeping only the first two or three 
terms respectively inside the curly bracket. We recall that exact integration 
of the truncated four (or three)-loop /3-function leads to a more involved 
structure than eq. (j3Tjj) . which poses serious difficulties in finding its inverse; 
nevertheless, in |I7j a useful solution has been still worked out at three-loop 
level via the real branch W_i(() of the Lambert function together with the 
Pade' Approximant of the related /^-function. Moreover, in ^J] a reliable 
approximation to higher order running coupling has been suggested, via a 
power expansion in the two-loop exact coupling eq. (j23J) . of the form 



V 



£ 

n>l 



# [a£V)]' 



(32) 



(k) 

with k > 3 the loop order, and p n proper functions of the coefficients (3j. 
Comparison of these multi-loop approximants to the coupling, with the 
relative asymptotic formulas, eq. (}3T]) and the three- loop analogue, reveals 
the better agreement of the former ones with the higher-loop exact coupling 
numerically estimated (i.e. starting from the exact implicit solution), even 
at low scales (see also |119j ). 

Finally, being the definition of the scaling constant A not completely unam- 
biguous, few comments are in order. As yet pointed out, starting from the 
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two-loop level an arbitrary constant has to be fixed for A being uniquely de- 
fined; beside the commonly accepted convention (|2U|) (or (|27|) ). we just men- 
tion that other prescriptions have been proposed 3U| in order to optimize the 
1 /L-expansion for the running coupling, while (}2*U|) does remain the preferred 
one as no further terms of order 1/L 2 appear in the two-loop asymptotic for- 
mula (J29j) . Thus, in the higher- loop levels the scaling constant is analogously 
related to the initial parameterization, and the four (and three) loop value 
reproducing the world average, roug hly a{M 2 z ) = 0.119, is A^" 5) = 220 
MeV with consistently five active flavors 0. 

A last remark concerns the scheme- dependence of the coupling and the 
scale parameter. Restricting ourselves to mass-independent RS (as MS-like 
schemes or trivially any prescription in the massless theory), renormalized 
coupling constants (see eq.((TJ)) in two such different schemes can be related 
perturbatively at any fixed scale \x 



a s /cts\ 2 
l + + v 2 — ) + 

7T V 7T / 



(33) 



Then, it is easy to verify that the first two coefficients of the relative (3- 
functions do not change as the renormalization prescription is changed, while, 
for instance, the third ones are related by (3' 2 = /3 2 —V\fi\ + v 2 j3o — (3qv\ . As a 
result the running coupling at each loop- level (e.g. eq.([3*Tj)) depends on the 
RS, through the coefficients j3j with j > 2 and the initial parameterization as 
well. The latter obviously amounts in suitably adjusting the scaling constant 
A , and the relation is given exactly by one-loop calculation [20] 

In (A'/A) = -\ , (34) 

which works through all orders; e.g. Af^g/A MS = exp [(ln47r — 7 e )/2] ~ 2.66, 
and with nf = 4 we roughly have [20] Amom/Ams — 4.76 where Amom 
refers to a scheme in which the subtraction of the relevant green functions is 
performed in the symmetric point and the renormalized coupling is defined 
through the three gluon vertex. 



2.3 Threshold matching 

Quark mass effects, till now ignored, reveal themselves in explicit corrections 
within higher order perturbation theory (see e.g. [21])) and in the energy 
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dependence of the effective (running) quark masses as a result of the RG 
improvement [22|, but we will not go into further details being these topics 
beyond the scope of this report. A direct effect of the quark masses on the 
evolution of the coupling is, as yet noted, through the dependence of the j3 
coefficients on the number rif of active quarks. A quark is active if rrif <C ix , 
where ix is the renormalization scale, and rrif the MS quark mass (see e.g. 
[18J); the definition can be also formulated in terms of the pole mass Mf that 
can be related to the former [23]. Indeed, within MS-like RS, decoupling of 
heavy quarks [21] is made explicit by constructing beside the full rif flavors 
QCD an (n/ — 1) effective theory below a heavy quark threshold [2H]. Then, 
to have a unique theory on the whole range, the two couplings ai n ^ and 
as ^ must be matched at each heavy quark mass scale li^^ — 0(rrif) . As 
a result the scaling constant A depends also on n/ (see e.g. [2H1 123 [28 j) . 
The most straightforward way is to impose continuity of a s by means of the 



matching condition at f ^(m 2 



fJ 



a 



m 



), which works up to next-to- 



leading order. At the one- loop level e.g., it translates into 



A (nf) = A (n/_1) 



A (n f -l)l 2/(33-2*,) 



rrif 



(35) 



which can be extended up to two-loop, and exhibits explicit dependence on 
the rrif values. Since trivial matching does not generally hold in higher orders 
within MS scheme, a more accurate formula is required in this case (see [21] 
and refs. therein); specifically to obtain the global evolution of the four- loop 
coupling the proper matching condition reads \27\ 



^( n / _1 ) A n f) 



l + k 2 



a 



(n f Y 



71 



a 



(njY 



71 



(36) 



with 



11 
72 



A;, 



564731 82043 



124416 27648 



2633 
31104 



(37) 



if /jW = rrif is exactly assumed. With this convention, eq.([3l)]) yields the 
relationship between the scaling constants A*™^ _1 ) and A^™^ in the MS scheme 
(see |27]). Note that one can equally fix /i*™^ = Mf , that amounts to a proper 
adjustment of the coefficients in (|3l)J): for instance [6J starting with (j3*T|) and 

A („ /=5 ) = 220 MeV, the values A^ n f=^ = 305 MeV and A( n /= 3 ) = 346 MeV 
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are obtained, with threshold fixed at the pole masses Mf, = 4.7 GeV and 
M c = 1.5 GeV. Finally, we observe that eq.(jHBJ) clearly spoils continuity of a s ; 
thereby, whenever continuity of the global coupling and of its first derivative 
is mandatory, one can resort to a more sophisticated technique [2H|, which 
relying upon mass-dependent RS yields a smooth transition across thresholds. 



2.4 Landau singularities 

The running coupling is not an observable by itself but plays the role of the 
expansion parameter for physical quantities. At one loop- level (cf. eq.(JH)), it 
resums an infinite series of leading-logs UV contributions and similarly the 
two-loop solution yields the next-to-leading-logs approximation; thus RG for- 
malism provides an iterative recipe for improving perturbative results. 
However, as a result, the RG resummation significantly modifies the analyt- 
ical structure of the series in the complex plane. As yet noted the one-loop 
coupling (fTHjl is clearly affected by a spacelike pole at A with residue l//3o 
(Landau ghost). Adding multi-loop corrections does not overcome the hur- 
dle. Rather the singularity structure of the higher order solutions is more 
involved due to the log-of-log dependence, so that a branch cut adds on to 
the one-loop single pole in the IR domain of the spacelike axis. Moreover, 
at a given loop level Landau singularities sensibly depend upon the approx- 
imation used. For instance, the two-loop iterative formula has a pole 
at z — 1 (/i = A , see eqs. fl^Tj) and with residue 1 / (2/3 ) , and a cut for 
< z < exp (— B\) due to the double logarithm. On the other hand, when 
considering the same loop approximation (J29|) . we note how the singularity 
in z = 1 already acquires a stronger character 

Po [z - I) 2 

with z = /i 2 /A 2 , and the cut now runs from to 1. Analogously the three 
and four- loop asymptotic solutions, as given by eq.(jSH), nearby the Landau 
ghost respectively becomes 

ro, N 5 2 m 2 (z-l) riv , £?? \n 3 (z — 1) , . 

^"i-trw- a - w "-Ai^TF (39) 



and are equally affected by an unphysical cut. However, the cumbersome 
singularity structure of the leading Landau ghost, and of the unphysical cut 
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as well, are an artifact of the UV approximations introduced. Therefore to 
deal with the IR behavior of the running coupling, e.g. at two-loop level, it 
is necessary to face with the exact solution (j23j): clearly it is singular when 
W_i(C) = -1 that is at z = B~ Bl (/i 2 = B± Bl A 2 ), corresponding to the 
branch point £ = — 1/e of the Lambert function. Nearby this point, due to 
the asymptotic (|25j) of W_i(C) , we have 



i.e., an integrable singularity (note that the awkward factor in front of A 
in the singular point can be reabsorbed into a proper redefinition of the 
integration constant, through (J2HJ) ) ■ 

A more detailed investigation about IR singularity structure of higher-order 
perturbative running coupling is performed on the ground of eq. ()32|) in the 
recent work |119j , where in particular the location of Landau singularities is 
determined as a function of rif . 

2.5 Time-like coupling 

Until now we have worked exclusively in the space-like region to derive the 
evolution of coupling, namely we have implicitly admitted the theory renor- 
malized at momentum scale with negative squared invariant mass. However, 
in perturbation theory one needs to parameterize observables depending on 
time-like arguments by means of an effective parameter, to improve per- 
turbative expansions. While this poses no special problem in large energy 
processes, at any finite energy the issue of which should be the most suitable 
parameter in the s-channel must be carefully considered. For this purpose 
we briefly sketch the key points of this subject, and start by noting that 
the standard practice is to merely take over to the time-like domain the 
space-like form at any loop level, regardless of the crossing between two 
disconnected regions, thus importing the same IR singular structure in a 
specular way. Nevertheless, from many early works based upon analysis of 
e + e~-annihilation data (see e.g. [3U] and refs. therein), it is known that this 
should not be the case but far in the asymptotic regime. This is because of 
the appearance of not negligible corrections (7r 2 -terms) to higher order coef- 
ficients of the a s (s)-expansions, due to analytic continuation from space-like 
to time-like axis. The problem has not yet been univocally solved, as it is 




(40) 
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strongly related to the IR non analyticity of perturbative running coupling, 
though it receives a well satisfactory answer in the framework of Analytic 
Perturbation Theory (see sees. 4. 3 and 4.6). 

Referring for defmiteness to e + e~-annihilation into hadrons, the issue can be 
stated as follows. Firstly recall that the main way (see e.g. [HI]) to deal with 
the ratio 

cr(e + e~ — > hadrons) . . 

= \ , _ i-zr- > 41 

cr(e + e — > fi + fi ) 

where s = q 2 > , is to start with the hadron contribution to photon polar- 
ization tensor in (space-like) momentum space 

nf (g) = OTg 3 - gV)lM-g 2 ) ex | d'xe^ < 0\T(f(x)f(0))\0 > (42) 

j 7 * being the quark electromagnetic current operator. Then, as known, op- 
tical theorem ensures the ratio R(s) to be straightforwardly related to the 
absorptive part of the forward scattering amplitude e + e~ — ► e + e~ . Indeed, 
being the analytical properties of the two-point correlation function (J4*2*|) in 
the cut complex plane C — {q 2 > 0} well established, this amounts to taking 
the discontinuity of Uh(—q 2 ) across the cut 

R(s) = lim [n h (-s + is) - U h (-s - is)] , (43) 
Ziyi e^o 

having computed the RG improved expansion for U h (—q 2 ) on the space-like 
axis (q 2 < 0). To this end one formally works with its first logarithmic 
derivative (thus avoiding subtraction constants), the Adler Z)-function [3*2*] . 
with the same space-like argument 



D( <r) >r- — • MO 

which in the MS reads 



,dU h (-q 2 " 
dq 2 



D PT (g 2 ) = 3^Q 2 



/ 

/ 



IX \ 7T / \ 7T 



being Q 2 = —q 2 , Q / the quark charges, and 
d 2 ~ 1.986-0.115n/ 



(45) 



d 3 ~ 18.244 -4.21671/ + 0.086nJ - 1.24 3 J^Qj ' ( 46 ) 
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Then by integration Uh(Q 2 ) is readily obtained. One should be aware that 
in the massless theory the cut spreads over the whole positive axis, and when 
taking into account quark masses the cut starts at the two-pion threshold 
4m?. . 

Whatever the loop order, the result for R(s) is usually recasted as a series 
in the effective parameter a s (s), naively obtained by specular reflection, that 
is by replacing the space-like argument Q 2 = —q 2 > straightway with 
the time-like one s = q 2 > in the coupling at a given loop order (e.g. 
eq. ()31Jl ). This final step displays nontrivial correction terms starting from 
0(«g), which are proportional to powers of it and are the drawback of the 
analytic continuation of hadronic tensor nearby the time-like axis. We then 
have the ordinary perturbative expansion for R(s) (see e.g. |Ho] ) 



R PT (s)=3j2Q} 



7T V 7T / \ TT 



(47) 



ir 2 b 2 

r 2 = d 2 ; r 3 = d 3 -5 3 , 5 3 = (48) 

with d 2 and d 3 the same as in 1)450. and we have used the shortcut bj = 
(47r) J " + /3j to emphasize the 7r-powers. The number 5 3 gives to the 0(a^) 
coefficient a strongly negative contribution for each n/, e.g. roughly ~ 14.3 
for rif — 4. Higher order 7r 2 -terms were analyzed in detail in j3HJ, for instance 
the fourth order correction is 

** = d*-U = ^(r> + gt-) (49) 



16 V 246 ^ 

roughly 5^ ~ 120 for n/ = 4; analogous (but more cumbersome) behavior is 
found for still higher orders, from which it becomes patent the remarkable 
growth in these correction terms. 

A similar treatment also holds for other s-channel observables as the 
normalized rate for r-decay into hadrons, showing that the effects of analyt- 
ical continuation make the perturbative expansions in the time-like region 
deeply different from Euclidean ones. 

Since the 7r 2 -terms play a dominant role in higher order coefficients, expan- 
sion l)17jl works only asymptotically at large s (that is when the smallness 
of a s (s) scales down these large coefficients); thus the space-like coupling is 
not a reliable expansion parameter in the s-channel at finite energy, and it 
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is not yet clear which one is instead to be used. Actually, as yet noted in 
pioneer works [HUEE], the expression of R(s), resulting from eq.(jIHJ) together 
with the improved perturbative series for 11^, exhibits no natural expansion 
parameter (since both real and imaginary parts of a s (Q 2 ) enter into the form 
of R), and the choice of such a suitable parameter is essentially a matter of 
expediency, that is it should be selected the one which yields better conver- 
gence properties. 

Alongside less meaningful attempts, we mention here the analysis [37| of the 
use of |a s (— s)| as expansion parameter for R(s) . By sensibly reducing higher 
order terms it yields faster convergence than a s (s); it further possesses the 
relevant feature of IR freezing, in agreement with contemporary phenomeno- 
logical models (e.g. [HH]), thus avoiding the hurdle of Landau singularities 
on the time-like domain. Indeed for low s , |a s (— s) \ is anyway less than 0.33 
for rif = 3 . According to this prescription the one- loop running coupling in 
this region should read [3*7j 



a s is) 



1 

Wo 



1 1/2 



ln 2 (s/A 2 



7T Z 



and asymptotically, i.e. for s ^> A 2 e" 



A)hi(s/A 2 



7T 



2 ln 2 (s/A 2 



+ 



(50) 



(51) 



resembling the UV behavior of the relative space-like coupling (fT%|) . Aside 
from these nice features, this function cannot entirely sum up the 7r 2 -terms. 
In order to deal with these corrections a somewhat different approach, known 
as RKP (Radyushkin-Krasnikov-Pivovarov) procedure, has been suggested 
jlOJ (see also [12]); which is firmly based upon the analytical proper- 
ties of the polarization tensor Il/ l (— q 2 ) and of the related D-function (J44j) . 
summarized by the dispersion relations, respectively 

R(s) 



n. 



ds- 



D{ 



ds- 



-q z 

R(s) 



v s-q 



2\2 



(52) 



(53) 



where R(s) is given by (J43|) . and q 2 lying in C — {q 2 = s > 0} . The key point 
here is the inverse of f!53j) given by the contour integral 

~ i£ dq 2 



ZlT 



(54) 
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to be computed along a path in the analyticity region for the D-function. 
Eq. (|54jl can be then generalized to an integral transform mapping a space-like 
argument function into a time-like one 



R(s) = $ [D(-q< 



(55) 



that can be straightforwardly applied to the perturbative expansion (|45p. 
provided that the integration contour is always kept far enough from IR 
space-like singularities. This yields R(s) as an expansion into the images of 
a s (Q 2 ) and of its powers, through the map <3> (10] 



R(s) 




n>l 



« S (Q 2 



7T 



(56) 



where d n are the same as in (|45|). Eq. ()56|) is to be compared with the stan- 
dard perturbative expansion ()47p: here the 7r 2 -terms are entirely summed 
up, with the drawback that within this framework there is no uniquely de- 
fined expansion parameter. However, it is useful to work out its behavior at 
0(a s )-approximation for R(s), i.e. 



$ [ai 1] (Q 2 )] =-<^-- -arctan 

fj IT 



1 1 



ln(s/A 2 



7T 



(57) 



easily obtained by applying the integral transformation (|54|) to the one-loop 
space-like coupling ([18)1 . The related time-like form (JSTj) is once more free of 
unphysical singularities at low s , and for s ^> A 2 e 7r as above can be expanded 
into powers of 7r/ ln(s/A 2 ) 



a (1) (s) 



7lp 



71 



7i" 



ln(s/A 2 ) 31n 3 (s/A 2 ) 



+ 



(5f 



Then eq. (l58|) can be recasted as a power series in the one-loop space-like 
coupling 



a s W( S ) 



1 - 



7r 2 6 2 /a^(s) 



48 



7T 



(59) 



emphasizing that the two couplings differ in three-loop level. By comparing 
(|59|) with eqs.(JHJ) and (|4*8|l . it becomes clear how RKP resummation of the 
7r 2 -terms into the time-like coupling (|57j) does work. The main shortcoming 
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of this recipe is that by applying to (JH7|) the inverse transformation to (|55jl . 
i.e. relation (JEHJ), we are not back to the original input (|18jl . Obviously, this 
is because integral transformations (p>4*j) and its inverse are well behaved as 
long as the integrand possesses the correct analytical properties in the cut 
complex plane; actually this is not the case for the space-like coupling (fTHj) 
and its higher-loop approximations, and we are then forced to compute the 
integral along a path large enough to avoid the IR space-like singularities. 



3 Infrared behavior 

As we told, Landau singularities are an artifact of the formalism and are 
related to a breakdown of the perturbative expansion even as an asymptotic 
expansion. On the basis of general principles of field theory a s (/j, 2 ) should 
have singularities in the fi 2 complex plane only on the negative real axis, from 
a threshold to — oo. This circumstance is irrelevant as long as we consider 
high energy processes in which the energy scale is much larger than A , but 
leads to serious difficulties in other problems in which much smaller values 
of \i are involved. Of this type are the bound state calculations, the decay of 
particles, the annihilation processes and even the deep inelastic collisions in 
particular geometries. Let us consider e.g. the ratio 

r(r ->• v T + hadrons) 
r(r -> v T + e + z/ e ) v ; 

between the hadronic and leptonic decay width of the r lepton. The non 
strange contribution to this quantity can be written in the form 

m 2 T J m 2 V m rJ 

l + ^)lmn2(t)+Imn2(t)} , (61) 

where Sew is the electro-weak factor, V u d the relevant CKM matrix ele- 
ment, n^j(t) and n^(i) are the transversal and the longitudinal part of 
the appropriate current-current correlator respectively. In this expression 
a s (/i 2 ) should be known in principle from the threshold ~ 0.14 GeV to 
m T = 1.78 GeV. As we told, to extrapolate a s (/i 2 ) to the infrared region, 
various proposals have been advanced. In this section we shall consider some 
purely phenomenological attempts, few examples of the so-called physical 
couplings and an optimization procedure. 
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3.1 Potential inspired approaches 

One of the first attempts to modify the expression for a s (/i 2 ) in the infrared 
region was made in the framework of the quark-antiquark potential. The 
potential V(r) between two infinitely heavy quark and antiquark (static po- 
tential) can be defined by the equation 

V(r) = Urn ^lnW[T] (62) 
T being a rectangular Wilson loop of size r x T and 

W[Y] = (~Tr7?exp Ug I dx»A^x)\ ) , (63) 



V being the ordering prescription on the path T and Tr the trace over the 
color indices. The simplest assumption is to write InW as the sum of its 
perturbative expression and a non perturbative term proportional to the 
area S = rT delimited by T 

i hi W = (i In W) VT + aS . (64) 

If (HnW)pT is evaluated at the first order in the coupling one obtains 

4 a s (/i 2 ) 

v { r ) = ~7T 1 - ar ; (65) 

3 r 

where /i should be taken of the order of the typical (1/r) . This is the so 
called Cornell potential, expressed as the sum of a Coulombian and linear 
term. As well known, by solving the corresponding Schroedinger equation 
a reasonable reproduction of the spin averaged bottonium and charmonium 
spectrum can already be obtained. In the momentum representation (|63|) 
can be written 

V(Q) = (k'\V(r)\k} = j^Jd 3 re-^V(r) = 



1 4a s (/i 2 ) 1 



a 



2vr 2 3 Q 2 tt 2 Q 4 



(66) 



where k, k' denote the initial and final center of mass momentum of the quark 
and Q = k' — k the momentum transfer 4 . According to the general rule we 

4 More properly stays for the limit of TQW^rp ~ jaz^rxs {the Fourier transform of 
—ir 2 re~ er ) for infinitesimal e. The same recipe should be applied to regularize eas. (|67l69(l . 
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should identify in (jSSj) the scale fi with Q. Since however in heavy quarkonia 
(bb and cc) Q typically ranges between 200 MeV and 2 GeV, we come close 
to Landau singularities and need to use some kind of regularization. In this 
order of ideas many years ago it was proposed [Hi] to write the entire potential 

as 

f/fn , 1 WQ 2 ) 

and to take 

«v(Q 2 ) = - ; .] Q ^ ■ (68) 
/3 ln(l + ^) 

Note that, defined in this way, ay(Q 2 ) reproduces the ordinary one-loop 
expression (JTSj) for Q —>■ oo, while for Q —>■ it gives 

A 2 

Replacing (p^l in (jfiTjl and comparing it with (jp^j) we see that (|5H|l incorporate 
the confining part of the potential (|65|) with 



«v(Q 2 )-^- (69) 



2A 2 , x 

(70) 

In (165)1 the perturbative part of the potential was evaluated at the first or- 
der in a s (/i 2 ). However, higher orders corrections have also been considered 
starting from (|6^jl or by other methods ^H]. As proposed in ^3], such cor- 
rections can be taken into account assuming (|67jl as the definition of a new 
coupling constant ay([i 2 ) and re-expressing this in term the ordinary MS 
constant. We can write 5 



a v (/i 2 ) = «is(^ 2 ) 

with 



TT \ TT J 



(71) 



31 

ui = — - —n/ ; w 2 = 28.538 -4.145^ + 0.077^ (72) 

i. e. u\ = 1.472, «2 — 13.190 for rif = 4. Alternatively one can use ordinary 
one or two- loop equations for the running ay(/i 2 ) but with A v redefined as 

(cf. « 

A v = A M ge^o ^lA2A m (73) 



3 Note that at three loops even terms of the form (^ a ) 4 lna^jg occur |46| . 
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and this should be the value to be used in Eq.(jHEJ can be re-obtained 
from (jl7j) by the modified /3-function 



fa (a) = -fa a 2 (1 - e^j . (74) 
More generally one can take 

"T^ + ^T^, (75) 



1 1 , Pi 1 .-la 



Z being an adjustable parameter. Eq. (|75|) reduces to (ff4^) for / — > oo, while 
for finite I produce the following asymptotic behavior 

fa(a) f3 a 2 - Aa 3 - . . . for a^O (76) 

fa(a) ~ — O! + . . . for a — > oo (77) 

which corresponds for ay(Q 2 ) to a two- loop expression of the type (|2b|) (or 
(jUD) for Q -> oo or of the form for Q -> . 

Eqs. (j67H7Uj) . or the more general corresponding to (|75|l . are very appeal- 
ing. They do not introduce any additional parameters in the theory and for 
Ay ~ 0.5 GeV (and I = 24) reproduce well the spin averaged first excited 
states in the bb and cc systems. However, they are not satisfactory from 
other points of view. They correspond to assume that all forces including 
confinement are due to the exchange of some effective vectorial object. On 
the contrary, if eq.([64p is generalized to Wilson loops with distorted con- 
tour (more sophisticate ansatzs also exist [SHI EI]) and some more elaborate 
method is applied, spin dependent and velocity dependent (relativistic) cor- 
rections can be obtained, which differ from those that would be derived from 
the exchange of such a vectorial particle alone and seem to be preferred by 
the data [17J 133 HH]- The majority of the analysis seems rather in favor of 
assuming a finite limit for a s (/x 2 ) for fi — > and adding a separate term for 
confinement as in (J65j) . Besides hadron spectroscopy there is other abundant 
phenomenology that seems to be consistent with the existence of a finite IR 
limit for a s (/i 2 ). This concerns the hadron-hadron scattering and the hadron 
form factors, the properties of jets, the transversal momentum spectrum in 
W and Z production etc. 

Coming to the explicit form of a s (Q 2 ) the simplest modification of the one- 
loop expression that has been considered is the "hard freeze" assumption 



for Q 2 > Q 2 



H = Wo MQg/A 2 ) tOT V < Wo 
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This equation has been used in hadron spectrum calculations, in a model 
for hadron-hadron scattering and in studies on nucleon structure function 
[!""""] (see also [SHI); the values Q = 0.44 GeV and A corresponding to 
— = 0.28 GeV have been found appropriate for the last two applications, a 
little smaller value - = 0.26 GeV comes from other phenomenology. Other 
convenient interpolation formulas between the large Q perturbative expres- 
sion and a finite a s (0) have been used again in hadron spectrum studies |54, 
with ~ 0.19 — 0.25. In a fully relativistic treatment in [S3] it was found 
that in order to obtain a ir mass so much lighter than the p mass a value 

= 0.265 was necessary. A similar results was obtained in j!65j with a 
one-loop analytic coupling (see later) with A = 0.18 GeV corresponding to 

= 0.44 but, what is really relevant, to an average value of - in the 
interval between and 0.5 GeV about 0.22. Finally in many analysis the 
successful phenomenological hypothesis was adopted that the gluon acquires 
an effective dynamical mass m g |57| 15**"] 139*] . We shall come back in the fol- 
lowing to this point, for the moment let us observe that to the leading order 
the following equation generalizes naively ("6*8*]) 



<Q 2 ) = — ^wv- ( 79 ) 

/3 ln(^^-) 

In particular the mentioned hadron-hadron scattering model gives in this 
case A = 0.3 GeV and m g = 0.37 GeV (^M = 0.26) jSH]. 
All the above mentioned attempts are essentially in agreement for what con- 
cerns the qualitative behavior of a s (Q 2 ) in the infrared region, even if they 
differ in the details. The point is that for many purposes we can simply 
introduce an infrared-ultraviolet matching point in the range of variability 
of Q; let us say Q\ = 1 or 2 GeV. Then one can use ordinary perturbative 
expressions of the type ("T~"]) for Q > Q\ and treat the quantities 

(80) 

or corresponding higher moments as adjustable parameters. Various event 
shape in e + e~ annihilation can be reproduced simply assuming Qi = 1 GeV 
and ( 2a ) ~ 0.2 [H3] . Such procedure is believed to be a way to parametrize 
and derive from experience truly non-perturbative effects (sec 4.5). 
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3.2 Physical couplings 

By physical couplings we mean any type of effective charges |U IHO] , defined 
by and in general any other observable related coupling. The advantage 
of such quantities is that they are in principle well defined at every scale and 
that can be easily extracted from experience. The disadvantage is that they 
are observable dependent. However, as we told they can be expanded in 
terms of the MS coupling (cf. |T2|) ) and related to each other [61"] . 
First we may consider the typical case of the quantity R e + e -(s) , (eq.(jJTJ)) 
and the related Adler function eq. (|44j) . Their perturbative expansions in the 
MS scheme are given by ()47|) and ()45|) . respectively. Correspondingly two 
different effective charges a^(s) and «d(Q 2 ) can be defined 

Dm=3j2Q} \i + a ^ 2 ' 

f 



/ 



Q) 



1 + 



a R (s) 



71 



TT 



The first of these is of time-like, the second of space-like type. As appar- 
ent from eq. (j53|l . they are related by the transformation (j55|l . i.e. = 
$ [«d(<5 2 )] , and are actually the model over which (JKTj) is written. 
A third effective charge, that has been used in analysis of the r decay, can 
be defined with reference to the quantity R T fcf. f!60|) and (jHU)); which we 
discuss in some more detail as an example. We can consider separately the 
contributions due to the vector and the axial currents and set 



127rS , Ew|Kd 

2t 

s 



dt 

s 



1 



1 



Imn^(t)+ImII$ A (t) 



L V/A* 



(82) 



For s = m 2 T eq. (|H2~|) gives the ordinary RX^ A for the r lepton, for s < m\ 



gives the corresponding expression for a fictitious r' with mass m T i = 
analogy with (|%T|) we can define a vector and axial coupling a^ A (s 



R^(s) 



2 



1 + 



V/A, 



s) 



71 



y/s. In 
by 

(83) 



where denotes the same quantity at zero order in the strong coupling. It 
is also defined a global a T (s) by 



R t (s)=RV(s) + rHs) = R < 1 



1 + 



7T 
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and a T (s) = \ [a^(s) + a^(s)] . The quantity Imlly'(s) may be assumed to 
vanish for small quark masses and ImlT^^s) to be given only by the pion 
pole, Imnl°'(s) = — ml). Under the same hypothesis the quantities 

Imily^(s) and Imll^(s) should be identical at the perturbative level and 
they are expected to differ asymptotically only for powers of 1/s which have 
a non-perturbative origin. The same must be true for (xf(s) and a^(s). 
Note that, due to isospin invariance Imlly'' is proportional to the isovectorial 
component Imllim 1 '' , and <xf(s) and «r(s) are related by 

Since the first coefficients of the perturbative expansion of R e + e - and R T in 
the MS scheme are known, even the coefficients Vi,v 2 ,v 3 of (j22J) and /3J, (31 
can be calculated. Therefore the couplings a T (s) and q;r(s) that can be ex- 
tracted directly from experience can be immediately translated in terms of 
a a^(s) value. Note, however, that in the expressions of t> 3 and /5 3 , which 
we do not report, there appears a quantity K± that has been only estimated, 
K± = 25 ± 50 (HE]- We report in fig|2fa) taken from the experimental 
value for a^(s), a^(s) and a T (s) as extracted from the data of OPAL col- 
laboration [66J. The results are confronted with the resolution of the RG 
equation for the approppriate (3 T (a T ) . As it can be seen the 3-loop a T (s) 
and 4-loop for K4 = ±25 have a finite limit for s — > ; the latter fits data 
very well for K4 = 25 down to s ~ 1 GeV . The strong enhancement of the 
experimental a T (s) below such value can be related to the pion pole that has 
not been included in the RG treatment. 

A last definition of effective charge has been given recently [69 in the con- 
text of the Sudakov resummation formalism. Let us consider, e.g., the Mellin 
transform F 2 (Q 2 , N) of the structure function F 2 (Q 2 , x) in DIS (x being the 
Bjorken variable). Even after the mutual cancellation of the infrared singu- 
larities due to the soft real and virtual gluons, F 2 (Q 2 ,x) has a logarithmic 
singularity at any order for x — » 1 which makes ordinary perturbation theory 
inapplicable. Such a singularity which is translated in the Mellin N-space in 
lniV power behavior for iV — > 00, can however be resummed [7U]. One can 
then obtain the following asymptotic equation 

d\nF 2 {Q\N) 

Q dQ 2 " ( } 
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9l 

71 



dx 



X 



N-l 



-A s [(l - x)Q 2 } + H(Q 2 ) + 



Uo 



where A$ and H have the usual expansions 

A S [Q 2 ] = a s (Q 2 )[l + a ias (Q 2 ) + ...}. (87) 

Then ^[Q 2 ] is assumed as an effective charge (Sudakov charge). For large 
rif the Borel transform of (|87|) -BfA^]^) can be given as an analytic expres- 
sion without singularities. The corresponding A^fQ 2 ] turns out to be free 
of Landau singularities but for Q — > behaves as — <^[^t — This em- 
barrassing singularity can be compensated with a non-perturbative term and 
leads to corrections in 1/Q (see sec 4.5). 

Another very interesting physical coupling is introduced in [BI]. This is a 
generalization to QCD of the Gell-Mann Low effective coupling for QED. 
It can be considered a generalization of the coupling ay(Q 2 ) discussed in 
sec. 3.1 and it is also related to the pinch scheme discussed in [71], which is 
somewhat more involved but explicitly gauge-independent. 
Let us isolate in a Feynman integrand the factor corresponding to the ex- 
change of one dressed gluon between two quark lines 

~^s(/i 2 ) 1 , x 

q 2 + t0 l-n[ g 2 ;a s (/i 2 ),/i 2 ]' 1 ) 



and set, e.g. in the MS scheme 



«sgd(Q 2 



Zt^u 2 ) 



qqg ' 



Zfq^iQ 2 ) 



l-n[-Q 2 ;a s (^ 2 )^ 2 ] 



(89) 



where Z^(fi 2 ) is the appropriate renormalization factor that makes the def- 
inition independent of the initial scale fi 2 . Obviously 



0^ a ) = [z a tf)z g (j?)]i = z m (v 2 )z;\v 2 ) 



(90) 



where we have used the same notation as in sec. 2.1. Note that in an abelian 
theory, like QED, Z^ = 1 due to Ward identity. In QCD the quark factor Z q 
cancels only the "abelian part" of the vertex factor Zq qg . This is the meaning 
of the superscript NA (non abelian part). The index SGD in eq.(|89|) means 
"single gluon dressing". 

From (|89|h if II is evaluated in the MS scheme, we can write, setting fi = Q, 



«sgd(Q 2 



a 



MS 



(Q 2 



a- 



MS 



(Q 2 



+ 



7T 



(91) 
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Here k\ depends on the constant terms occurring in the renormalized II and 
it is so gauge dependent; e. g. in the Feynman gauge k\ = (31 — 10n//3)/12. 
Note, however, that eq. ()89|) can be related to the pinch scheme and k± 
made gauge independent simply including in the definition of II additional 
constants coming from the pinch parts (selfenergy like parts originating from 
contraction of internal lines as a concequence of Ward identities) of the vertex 
and the box diagrams [71]. In this way we obtain k\ = (67 — 10rij-/3)/12 . On 
the contrary, if we absorb the pertinent soft gluon corrections in a redefinition 
of a S GD(Q 2 ), we find h = (67 - 3vr 2 - 10n//3)/12 [221 EH- It is the latter 
value which is used in connection with the considerations of sec. 4.5. 



3.3 Optimized perturbation theory 

By optimized perturbation theory it is generally meant some kind of per- 
turbative expansion in which the expansion variable, or the splitting of the 
Lagrangian in an unperturbed and a perturbation part, is chosen in depen- 
dence on a number of arbitrary parameters. Such parameters should not 
appear in the exact result, but obviously they occur in any expansion which 
is stopped to a certain maximum term n. However, it is immediately shown 
even by very simple examples, that the convergence of the series is greatly 
improved if for every n the parameters are chosen at some stationary value 
(optimized choice) that depends on n. 

In QCD, due to the arbitrariness in the choice of the RS, the optimization 
can be required in the variables that control such a scheme, e.g. the subtrac- 
tion point fi and the scheme dependent coefficients /3 2 , /?3, • • • • Obviously 
every choice corresponds to a different definition of the coupling constant a s . 
Let us e.g. consider the quantity R e + e -(s) defined by eq. (}4~T]) . and rewrite 
expansion (|47J) with an arbitrary choice of /i 2 (for the moment different from 
the total energy s) and in a arbitrary RS 



R e +e- 



2 
f 

f 



i + *w + ,, 2(s) ( <>m 2 + rsW f<tLi\ +... 



7T \ 7T / \ 7T 



(92) 



The quantity R e + e - (s) must be RS independent and, if we neglect the masses 
of active quarks, we can write 
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d f as a' j+2 d , 

P{a s ) / da'j—-^— R e+e - = , (93) 



where j — 2, 3, . . ., r = ln(/i 2 /A 2 ) (recall eq.(|2%J)) and we have used 

da 1 __ R( s f aB a ' a " J+2 - Q s +1 ( 1 A J" 2 

(94) 

as it can be seen deriving eq.(JT7J). Eq.(|9*3^l can be used, first to obtain 
r 2 , t 3 , . . . in an arbitrary RS when we know this quantities in a specific 
RS (cf. eq.flnj)), and then to make the optimal choice for r, fa, fa, ■ ■ •, when 
we have stopped expansions (JH2*|) and (JHJ) to a certain maximum order; let us 
say to the term in r 3 in (}9*2"j) and to three- loop (i.e. to the term in fa ) in 
From now on we will use the notations bj = (47r) J+1 /3j and a = a s /ir . Re- 
placing (J92*Jl in (J9~3~j) . asking that this equations for a given a 2 are satisfied 
for an arbitrary value of a we obtain differential equations for r 2 , r 3 , . . . . 
Restricting to r 2 and r 3 and j = 2 we have 

dr 2 _ 1 <9r 2 _ 

W = W 2 + B bl 1^ = -^ (95) 
Integrating the above equations, we obtain 

T2 = ^b r + p 2 

1 l2 2 1, 1 , 1 h 

r 3 = — o n r H — OoP2 r H OiT — h p q 

3 16 2 0F 16 16 6 H3 

IbA 2 1 b 2 

r ' + si) (96) 

where p 2 and p 3 are integration constants and so quantities independent of 
r, b 2 , ... and RS independent. They can be calculated e.g. equating b 2 , r 2 , r 3 
to their expressions b^ 5 , r^ s , r^ s in the MS scheme as given by equations 
and (|46|) after setting u 2 = s . We have 
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Note that p 3 turns out to be independent of s , and r 2 has the form 



r 2 



— -b In ^— + -6 r + t 



A1S 



4 



A 2 4 



(91 



while r 3 depends on s and r only through r 2 . Let us now replace /3(a) in (|17|) 
with its 3-loop expression P^(a) 



-7ra 2 (| + fa + ffa 2 ). We have 



4 , 6i /6 x a 



&i , f 166 + 46ia + 6 2 a 2 
m 



26 2 



, 26260-6? 
+ — /(a, 62) 



26 2 



with 



/(a, 6 2 



1 . 46 
in ■ 



ia(6i + v^D) 



(99) 
(100) 



46 + \a{b x - y/D) 

and D = b\ —46260 . Let us make the same replacement in (J93)) and stop (J92|) 
at the a 3 term. By requiring that eq.([93p is now exactly satisfied, we obtain 
the following equations 

,2 



36 r 3 + ^6ir 2 + \-b 2 + (36ir 3 + -62^)7 + 36 2 r 3 ^- = 
I lb z 4 lb 



46 



2r, 



/(a, 6 2 ) — a = 



with 



1(0,62) = -jj- 



(46 2 - Sb 2 b )a + b 2 b 1 a 2 
lb6 + 46 x a + 6 2 a 2 



2b 2 b f(a,b 2 



(101) 



(102) 



Eq. (j9*9"|) gives a as a function of r (or /i) and then ()102|) and 1)101)) become 
equations in r and 6 2 that determine the optimal choice r(s) and 6 2 (s) of such 
quantities for every s. In this way we obtain an optimized running coupling 
which together with the optimized values f 2 (s), r 3 (s), can be 



as 



Qqpt(s) 



used in (|9*3*jl to evaluate the quantity i? e + e -(s). The resulting aopT(s) can be 
evaluated numerically and is reported in figJ2fb) taken from [72] as a function 
ofq = y/s for A m = 230 MeV (Ams = 264 MeV). Note that a OPT (q 2 ) stays 
finite as q decreases and attains a maximum for q ~200 MeV, after which 
it remains practically constant. Such a maximum value can be shown to be 
given by the equation [72] 
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Figure 2: (a) Effective charge a T compared with solution to RG equation up to four- loop, 
(b) The optimized third-order results for the coupling a = a s /ir and the QCD correction 
TZ^ to the ratio R e + e - ■ Also shown is the second-order result TZ^ . Quark thresholds are 
indicated by vertical lines. 



Taking n/ = 2, as appropriate to the range of energy, we have from (jHZJ) P3 = 
-10.911 and so = 0.263. Note that the behavior of a OPT (g 2 ) for < 
q < 1 GeV is perfectly consistent with what follows eq.()8L)jl. In figJ2fb) are 
also reported the second and third order results for the quantity 71 defined by 
R e + e - = 3 <5/(l + 7V) and expressing the QCD correction to R e + e - . The 
results compare favorably with the experimental data appropriately smeared 
to wash the irregularity due to the resonances. 

4 The dispersive approach 

As mentioned in the introduction we can dispose of the Landau singularities 
simply exploiting the general analyticity properties expected for a s (Q 2 ) and 
applying perturbation theory directly to the spectral function (7^1 IHDJ ■ This 
idea generalizes to QCD a method originally introduced in QED [HI]. 
Defining the spectral density 



p(cr) = Ima s (-<7 - iO) = — [a s {-a - iO) - a 8 (-a + i0)] , (104) 
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where a > and a s (— a) is the perturbative RG solution at a given loop 
level, the analytically improved running coupling is thus given by 

«an(g 2 ) = - r d(J -^k ( 105 ) 

whose argument Q 2 = —q 2 > now runs over the whole space-like axis (we 
observe here and in the foregoing the identification /x 2 = Q 2 ), that is a an (Q 2 ) 
is free of any space-like unphysical singularities by construction; moreover, 
due to the asymptotically free nature of the perturbative input the spectral 
integral needs no subtractions. 

Note that different strategies to incorporate analyticity into the RG formal- 
ism, or even to implement the above device, exist as well, and they will be 
briefly reminded in sec. 4.4. 



4.1 One-loop analytic coupling 

At one-loop level this trick works quite straightforwardly; starting from the 
leading-logs expression (fTH|) the related spectral density (see also eq.([13|)) 



\n 2 (a/A 2 ) + 7i 2 



(106) 



turns out to be explicitly integrable, and eq. ()105|) yields 



+ 



A 2 



_ln(Q 2 /A 2 ) A 2 -Q 



(107) 



The analytically generated non-perturbative contribution in ()107j) subtracts 
the pole in a minimal way, yielding a ghost-free behavior which avoids any 
adjustable parameter. Obviously eq.(jHJ) for the scaling constant does not 
work anymore and, at one-loop, it has now to be changed to 

A 2 = /i 2 exp[-0(/3 o asG"o))] , (108) 

where the function is related to the formal inverse of ()107)) that is, with 
x = /9o«an and Q 2 /A 2 = exp0(x), it satisfies 



4>(x) 1 — exp <p(x) 
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Among the main features of the analytically improved space-like coupling, it 
should be firstly stressed its agreement with asymptotic freedom constraint, 
being the pure perturbative contribution ruling in the deep UV region over 
the "non-perturbative" one. Indeed the latter for Q 2 > A 2 can be rewritten 
as the sum of the series 

1 00 /A 2 \ n 

■ (no) 

yielding power correction terms to the one-loop perturbative coupling |T8j) 
(see sec 4.5). On the other hand, in the extreme opposite domain eq. (jl()7)) 
exhibits the infrared freezing value ain(O) = 1/(3q — 1-396 choosing consis- 
tently nf = 3 at low scales; this value turns out to be independent of A and 
universal with respect to higher-loop corrections, i.e. the analytic coupling 
()105|) has a remarkably stable IR behavior. 
The beta-function for the one-loop coupling ()107|) reads (73] 



' kI {x) [exp <j)(x) - 1] 



with satisfying eq. (jl()9j) . Despite the implicit form of (jlllj) . its symmetry 
property /3^(x) = (3^(1 — x) reveals the existence of a IR fixed point at 
x = 1, corresponding to a au (0) = l/Po ( scc a l so [88J) . 

Finally, we just mention here that the relation between the /^-function struc- 
ture and the analytical properties of the running coupling has been inves- 
tigated in ref. [HI] within a pure perturbative approach; some resummation 
tricks for the asymptotic /^-function have been there outlined, which can 
cure the Landau ghost problem, and rely upon freedom in choosing the RS 
in higher orders. Another viewpoint is given in J7j, where the hypothesis of 
pure perturbative freezing has been explored by analyzing the nj-dependence 
of the /^-coefficients (see also eq.(|S7|) and the following). 



4.2 Two-loop and higher orders 

Actually, as discussed in sec. 2.2, the two and higher-loop RG equation for 
invariant coupling has no simple exact solution, so that rough UV approxima- 
tions are commonly used (eq.([31|)). leading to a cumbersome IR nonanalytical 
structure. Nevertheless, to go further inside analytization and its features we 
follow ref. |82| and choose the two- loop iterative solution as given by eq. (|26j) . 
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By applying now the analytization recipe the main difficulty arises from the 
integration of the related spectral density [ZHj 

^^=^^ = hmfm ^ Ha/A2) (112) 

I(t) = n + Bl arccos /( Bl + J \ , R(t) = t + B x In ^ + g + ^ 

y/(Bx + t) 2 + 7T 2 Si 

which does not lead to an explicit final formula, though the relative dispersion 
integral (jl05|) can be handled by numerical tools. However, recalling the 
singularity structure of (sec. 2.4), entirely subtracted by analytization, 
the two-loop analytic coupling can be recovered by merely adding to (J25J) 
two compensating terms, cancelling respectively the pole and the cut 



a®. t ( z ) = a ®(z) + AW + AW (113) 
A?>(*) 



p w 2/5 1 - 2 

1 fwl-*) ^ £ x 



A) Jo £ - 2 [lnf + Si In (-1 - Bf 1 lnf)] + tt 2 S 2 

with the dimensionless variable £ = a/A 2 . Despite the little handiness of 
(|113|) . one can still readily verify its limit ot^ At (0) = 1/A)j f° r more general 
arguments concerning universality of the IR freezing value through all orders 
see for instance |H2"1 IHB*] and [HZ1- Moreover, the non-perturbative UV tail 
of analytized coupling can be estimated by expanding the two compensating 
terms in (|113|) into inverse powers of z = Q 2 /A 2 for large z [85] 

A( 2 )(,) + A( 2 )(,) = l^| (114) 
1_ exp[-nB l( l + 0] 



2 A '(l + £-ln£)' 



^ 2 



Thus for z > 1 non-perturbative contributions can be recasted analogously to 
()110|) as a convergent power series with all negative coefficients, approaching 
-1/2 (e.g. ci = 0.535, see refs. [HS] and [EE]). 

As yet pointed out, to get the most accurate result at two- loop in the IR 
domain, one needs to start with the exact RG solution (|23|) . To this aim first 
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of all one has to extrapolate eq. (]23[) in such a way to define a real analytic 
function regular for — e _1 < £ < 0; the two branches Wi(Q and W_i(0 
merge continuously on (— e^ 1 , 0) , so that the correct recipe is to use eq.(j23J) 
for Im(C) > and the same equation with W_i(C) replaced by Wi(() for 
Im(£) < . The discontinuity of the two loop exact solution across the 
time- like cut defines the spectral density jHO] (see also [BE]) 



C(t) 



-Im 



eB 1 



cxp 



Li + wi(c(t))J 



h 11T [ 1 

B 1 \B X 



ln(<x/A 2 



;ii5) 



Dispersion integral with t = ln(cr/A 2 ) now leads to the exact two-loop 

analytic coupling 

i r° . e* 



a 



(2) 



7T 



dt 



e l + z 



P {2) 

rex 



(116) 



where pix (^) = pix (cr) . Numerical estimates of (fTTBj) as well as for the 
analytic iterative coupling ()113|) have been performed at low scales with 
n,f = 3 in [HE] (see Tab.l), both normalized at the r mass M T = 1.777 GeV, 
a s (M 2 ) = 0.36; comparison reveals the relative error for the analytized solu- 
tion ()113J) to be around 1.8% in the IR region. 

In spite of its accuracy eq. ()116|) cannot be easily handled, even though it is 
a source of numerical information to which simpler expressions have to be 
compared. Actually, for practical aims, many useful two-loop approximate 
formulas have been suggested; among them we recall here the "one-loop-like" 
model jE2 



a 



(2) 



(0 



1 

A) 



l 



exp / — 1 



I = In 



A 2 



B x In Win 2 



91 
A 2 



47T 5 



(117) 

suitable for analysis of rather low energy phenomena, since it approximates 
the exact two- loop analytic coupling with 1% precision for Q > 1 GeV, and 
correctly reproduces both the universal freezing value and the UV two-loop 
asymptotic behavior. However, its accuracy breaks down when taking into 
account flavor thresholds. To this end it has been suggested |117j to use 
eq. ()117|) provided that the scaling constant and the coefficient B\ are re- 
placed by adjustable parameters (respectively A 2 i and c^f listed in tab. Ill 
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of ref. |117j for different initial A^ n f =3 ^), as a result of an interpolation pro- 
cedure; this yields an accuracy within 1% in the whole space-like region. A 
"one-loop-like" model employing the scaling constant as a fitting parameter, 
suitable from 2 to 100 GeV, has been recently developed in [55] . 
Increasing difficulties arise when dealing with even higher-loop level and it 
becomes prohibitive to achieve useful explicit formulas. Starting e.g. with 
the standard three or four-loop asymptotic solution (|31j). one has to face 
with the leading singularity in z = 1 of the form (jHHj), beside the IR log-of- 
log generated cut; terms accounting for these divergences acquire the form 
of cumbersome finite limits integral as in ()113j) . Nonetheless, the effects of 
non-perturbative contributions have been widely investigated up to four-loop 
jSSllHZl, both in IR (where they play the most prominent role) and UV region, 
by using asymptotic solution (j3"Tj) as a perturbative input. While confirming 
at once IR stability due to the universal freezing value of the analytized cou- 
pling, its UV tail has been reduced jHZj in the form of power type corrections 
analogous to IJ114j) . Within this approximation the n/-dependent coefficients 
c n are all negative up to four- loop, and their absolute values monotonously 
increase with n. In the large Q 2 limit, however, there is no need to sum a 
high number of terms, and truncation of e.g. the three-loop non-perturbative 
series to first term yields the approximate expression [87J 

Po V 



2 



B 2 

ci = -l + Bi(l-7te)--^ 



- 2 - . 2 



B 2 - — + (l- 7e) 
o 



for the analytic coupling, with 1% accuracy yet at Q 2 ~ 5 A 2 . Here 7e is the 
Euler constant, B 2 = ftofa/ ftf and ai 3 \z) is the perturbative counterpart as 
given by (pTTj) ; if rif = 6 we roughly have c\ ~ — 0.52 . 

By normalizing the three-loop analytized coupling at a s (M|) ~ 0.118 (with 
perturbative input from eq.([31|) and numerical evaluation of dispersion in- 
tegral HEED ), one can extract the three-loop scaling constant A^ 5 ) ~ 
210 MeV jHZj, which lies within the errors of the pure perturbative esti- 
mate (see sec. 2.3), being the non-perturbative tail negligible around the 
normalization point. Obviously the main discrepancies emerge in the low 
energy region, where non-perturbative contributions slow down the rise of 
the curve. By using continuous matching up to three-loop, at the MS quark 
masses m;, = 4.3 GeV and m c = 1.3 GeV, for the three-loop analytized case 
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one finds roughly jHZj 

AW=5) ~ 210 MeV -> A (n ' =4) ~ 299 MeV -> A (n / =3) ~ 382 MeV, (119) 
to be compared with the relative perturbative estimates 6 

A (n ' =5) ~ 210 MeV -> A (nf =4) ~ 290 MeV -> A (n/ =3) ~ 329 MeV . (120) 

These discrepancies at low scales can be translated into the values of analytic 
and perturbative coupling at the r mass M T = 1.777 GeV, given respectively 
by oiX M r) - °- 294 and al 3) (M T 2 ) ~ 0.318. Note that normalization at 
a s (M T ) = 0.35 is sometimes adopted, and this leads to considerably higher 
values for the scaling constant. Last, it should be stressed the nearly indis- 
tinguishability of the two to four- loop analytic curves (figG^a)), confirming 
higher order stability on the whole spacelike axis [751 Ell- Thus we argue the 
two-loop correction to be a well satisfactory improvement of the one-loop 
result, so that it is reasonable to resort to approximate formulas as given e.g. 
by (|117|) . However, in recent works |9*T| I119j the multi-loop approximation 
(|32|) as a power series in the two-loop exact coupling (J23j) has been exploited 
as an input in the dispersion integral ([105)1 . to yield high-accuracy three and 
four-loop analytic coupling of the form |5T] 

«i fc n ) (Q 2 )=E^ ) K 2 n ) (^)] n (121) 
n>l 

with ain given by ([116)) and coefficients the same as in ()32[) . Numerical 
values at low scales can be found in [9"T] . to which we refer for more details. 

4.3 The time-like coupling in the analytic approach 

Coming back to the issue of defining a reasonable expansion parameter in 
the time-like domain (q 2 = s > 0), we start by noting that such a definition 
naturally arises in a self-consistent way within the framework of the analytic 
approach Further, it can be regarded as the final step in the proce- 

dure RKP for 7r 2 -resummation outlined in sec. 2.5. 

As yet noted, eq.(j53J) and its formal inverse (jo"lj) can be generalized to the 

6 The shifts of As in 1|120|) w.r.t. sec. 2.3 are mainly due to the different normalization 
value, and partly to continuous matching used here; indeed there is little sensitivity to the 
implementation to trivial matching given by ea. (136(1 . 



42 



proper tool for relating s- and t-channel observables, that is for crossing the 
two distinct regions. On this ground, it has then been suggested |12 E2j to 
use the same integral transformations in order to connect the time-like and 
space-like couplings, employing as a suitable regularization the analytically 
improved coupling ([105)1 . free of unphysical singularities at any loop level. 
Then |22| 

i f s+ie da 2 

«( S ) = W ^«an(-g 2 ) (122) 



2 

oo 



« an (-g 2 ) = ~q 2 I ds a[S) (123) 
Jo ( s - 1 ) 

with q 2 as usual in C — {q 2 = s > 0} . Specifically, eq. ()122|) can be assumed as 
a definition of the time-like coupling. Once the space-like singularities have 
been washed out the ambiguity about the path for ()122|) disappears, and we 
find a one-to-one relation between the t- and s-channel couplings. Note also 
that a(s) can be equally defined by the differential equation 

s —a(s) = --p(s) with <5(oo) = 0, (124) 

as 71 

as can be immediately checked, e.g. by differentiating eq. ()123|) . Then, 
eq. ()124j) immediately yields 



as 



1 f°° da , N . . 

/ —M (125) 



and p(a) as given by ()104|) . Moreover, eq. p24|) emphasizes the straightfor- 
ward relation between the "time-like /3-function" and the spectral density, 
thus reviving an old hypothesis due to Schwinger [Mj . 
With the use of ()124|) . eq. (|123|) can be also formally inverted as 



sin(7r'P) 1 2 

a(s) = — — a an (s) = a an (s) - — {irV) a an (s) + . . . (126) 

71 F 3! 



a an (s) 



7T 2 (3 2 2 57T 2 /?o/3i 3 

—z-a m (s) aL(s) + ... 



where V = s(d/ds). Of course at one-loop ()126|) coincides with ()5T?j) . and 
from eqs. ()125|) and (|106|) one finds again eq. ()57|) : however this now leads 
through (|123|) to the starting space-like coupling (|1(J7|) . being the required 
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Figure 3: (a) Analytic coupling from [57] up to four-loop with the above normalizing con- 
ditions and threshold matching, compared to pure perturbative input l|31l) . (b) From[53]: 
comparison of one-loop timelike and spacelike couplings i|57[l and (|107|l . 

analytic properties preserved within this framework (they are both plotted in 
figEKb)). Furthermore, at two- loop level inserting into eq. (|125|) the spectral 
density (|115J) computed on the two-loop exact RG solution (|2^|). unlike the 
space-like case, integral can be taken analytically 



COO 



Pi n(3 



Im[l+W! (C(s))] 



(127) 



eB l 



exp 



VK 



1 



— - 1 



The main feature [92] of ()125j) and ()105j) is the common freezing value at 
the origin a(0) = a an (0) = l/fio, independent of the loop level and of any 
adjustable parameter. Moreover, they exhibit similar leading UV behavior, 
constrained by asymptotic freedom, as can be seen by ()126j) (and recalling 
that asymptotically ()105|) reduces to the pure perturbative coupling a s ). 
Nevertheless this approximate "mirror symmetry" is broken in the interme- 
diate region, the discrepancy being about 9% at one-loop, and slightly less 
at two and three-loop (see ref.jH21 f° r numerical comparisons). 
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We finally just note that ref-jHE] exploits an argument against a possible ex- 
act symmetry ruling the "t-s dual" couplings, ()l()5j) and ()125j) . on the ground 
of causality principle. 

4.4 Related models 

Since the last ten years there have been a number of different attempts to 
avoid Landau singularities invoking as well analyticity of the coupling in 
the space-like momentum region. Within the dispersive approach, it is re- 
markable the existence of models suggesting IR enhancement of the QCD 
coupling, whose most attractive feature is supposed to be a straightforward 
relation with quark confining potential within the framework of one-gluon 
exchange (see sec. 3.1). 

To this class belongs for instance the "synthetic coupling" model recently 
developed in [HHl which amounts to modify by hand the analytically 
improved coupling (|1(J5|) by additional non-perturbative pole-type terms; at 
one-loop it reads [HZ] 



a syn (Q 2 ) = 7T 

Po 



A 2 cA 2 (1 - c)A 

+ 77 7^ + ^TT + 



ln(Q 2 /A 2 ) A 2 -Q 2 Q 2 Q 2 + m 



(128) 



where m g = A./y/c — 1 . Infrared enhancement, due to the pole term at 
Q 2 = (cf. eq.(p9^l). is governed by one dimensionless parameter c G (1, oo), 
which is meant to relate the scaling constant A to the string tension o of 
potential models, through the same eq.(J7UJ) but with A replaced by -Jch. The 
pole term at Q 2 = —m 2 < corresponds to a non vanishing dynamical gluon 
mass, while leaving the analytical structure of eq. (|107|) along the space-like 
axis unchanged. The value of the m g parameter has been estimated [HZ] as 
400 - 600MeV (see sec. 3.1). Ea. dT^ can be derived, analogously to (fTUTj) . 
from a dispersion relation with a spectral density of the form (J 106)) plus two 
5-terms properly accounting for the poles. Along with the singular behavior 
as 1/Q 2 at Q 2 — > , reproducing the linear confining part of the potential 
(jBT|) . construction of ()128)) is mainly motivated by the UV asymptotic jHZ] 
of its non-perturbative contribution of the form 1/(Q 2 ) 3 , decreasing faster 
than ([TTTIjl as Q 2 ->• oo . 

Quite analogous result has been achieved in [HB1 IHE] , merging analyticity and 
IR slavery at zero momentum in the one-loop formula 



(129) 
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where z = Q 2 /A 2 . The trick undertaken here is to impose analyticity ab 
initio on the whole perturbative expansion of the QCD /3-function, and then 
to solve the ensuing " analytized RG equation" for the running coupling. 
This could be done in principle at any loop-level, formally 



where the r.h.s. of ()130|) is achieved through the usual dispersion integral 
(|105|) . starting from its spectral density now given by the discontinuity of 
the expression of the 1-loop /3-function as a whole across the time-like cut. 
Note that in the one- loop case the r.h.s. of (|130|) is merely proportional to 
the one-loop analytic coupling f)407|) . so that further integration leads to the 
singularity at zero momentum. This "new analytic invariant charge" pos- 
sesses the universal (i.e. loop-independent) asymptotic ~ A 2 /Q 2 as Q 2 — > , 
which again results in confining quark-antiquark potential [§H|. Besides, this 
variant of the dispersive approach, while disclosing in some measure the am- 
biguities suffered by the method, shares appealing features with the IR finite 
analytic coupling (jl05|) ; namely, it displays no unphysical singularities and no 
adjustable parameters, and exhibits good higher-loop and RS stability (see 
e.g. fE^ for technical details). Analogous IR behavior of the QCD coupling 
has been found in [100] on the ground of different reasoning (see also [HHl EH] 
and refs. therein for similar proposal and results). It is worth noting that 
in the most recent developments |lUlj of the model in hand, eqs. ()130p and 
(j!29|) . inclusion of the lightest hadron masses (ir meson) is accomplished con- 
sistently with the dispersive approach, and relations with chiral symmetry 
breaking phenomena are also investigated |lU2j . It has been shown there how 
nonzero pion mass substantially affects low energy behavior of the invariant 
charge f)129j) . by slowing down the IR enhancement of the massless case to 
an IR finite value, depending in turn on the pion mass. As a result, this 
massive running coupling displays a plateau-like freezing on the IR time-like 
axis, specifically for ^fs in the interval between and the two-pion threshold 
(see also |l()3j ). in qualitative agreement with results of OPT (22j (sec. 3.3). 
Finally we remind another attempt |lU4j to modify eq. (jlU7|) in order to es- 
timate non-perturbative power corrections (see sec 4.5). This is dictated by 
the need to cancel somehow the unwanted behavior (jllOJ) . and further to slow 
down the too large value of (jl07|) at the origin (see sec. 5). This has been 
done by requiring no power corrections faster than 1 / Q 2p , and a number of 




(130) 
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adjustable parameters to remodel the IR behavior of the one-loop analytic 
coupling (jl()7j) . Thus an useful generalization reads 



«w0) = 7T 



I z + b fl + c 



\nz (1 - z)(l + b) \z + c 



in 



(131) 



Setting b = 1/4 and c = p = 4 and n f = 3, A = 250 MeV eq. lfHCT) has 
a maximum at 0.4 GeV and then freezes to a considerably lower value than 
([107)1 . fitting data on power corrections (see |104j and refs. therein). 

4.5 Power suppressed non-perturbative corrections 

As mentioned in the introduction, intrinsically non-perturbative effects would 
manifest themselves in power-type corrections A u /Q u or (A u /Q u ) \n(Q 2 /Q 2 ), 
to the expression of various observables. We want to discuss the subject in 
the context of ref. [HI], in which it was originally proposed, by making use 
of the coupling (Xsgd(Q 2 ) we have introduced by eq. ()89|). 
Note that, as consequence of its definition, o;sgd(Q 2 ) must have only physical 
singularities in the Q 2 complex plane and the entire dispersive formalism 
discussed in this section can be applied to it. Actually such formalism sprang 
from a common source |105j and developed along parallel lines from references 
[75] and jOl]. According to eq. (|105|) we can write 

= Q *» (0) -v/ + cpy (132) 

Then, the factor corresponding to a dressed gluon line in a Feynman integral 
can be written 

-^asGD(-^ 2 ) -««sgd(0) 1 f°° dm 2 2 -i 

"P( m )75 ( 133 ) 



fc 2 + zO A; 2 + iO 7i J m 2 k 2 — m 2 + iO 

Let us now consider some hard process initiated by a quark and an inclusive 
infrared and collinear safe observable V(Q 2 , x) related to it. Let us assume, 
to be specific, the observable has a zero order expression Vq(Q 2 ,x) in terms 
of parton model and consider the first order QCD correction Vi(Q 2 ,x) in 
which dressed gluon lines are inserted in the original skeleton graph. Let us 
denote by T\{e, x) the Feynman integral (or the sum of Feynman integrals) 
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that gives such correction but in which formally a mass m is given to the 
gluon and it is set e = §5-. Due to (|133|) the insertion gives 

1 f°° dm 2 

^(Q 2 ,x) = a S GD(0)^i(0, x)-- p {u?)T x {e, x) = (134) 



it j m? 



r°° d m 2 . . Q 

= / — 5-asGD(m 2 )J'i(e, x) , with T\{e, x) = -e— T\{e, x) , 
J m z oe 

where we have introduced the time-like coupling 7 «sGD( m2 ) , related to 
«sgd(^ 2 ) by the equation m 2 {da^Q^m 2 ) / dm 2 ) = —p(m 2 )/n (cf. eq. ()124|) ) . 
Eq. ()134j) gives a kind of weighted average of !Fi(e, x) and shows that the 
effect of the running coupling can be simulated giving an appropriate effec- 
tive mass m g , and justifies the phenomeno logical application based on such 
idea that we mentioned in sec. 3.1. Let us now write q;sgd(<5 2 ) an d corre- 
spondingly <Ssgd(<3 2 ) as the sum of a perturbative and a non perturbative 
part 

«sgd(Q 2 ) = a^ D (Q 2 ) + a^ D (Q 2 ) • (135) 

Since, however, perturbative theory appears to work very well down to Q ~ 
1 or 2 GeV and in some particular case even better, o^q D (Q 2 ) must vanish 
sufficiently fast as Q increases, let us say at least as l/Q 6 . Consequently, 
looking at eq. (|123|) we must assume that at least the two first integer moments 
of <Ssgd(^ 2 ) vanish, 

f ^m 2 ^ G p D (m 2 ) = 0, r ^m 4 ^ G P D (m 2 ) = 0. (136) 
Jo m Jo m 

Let us now restrict to collinear and infrared safe observables. Then !Fi(e, x) 
must vanish conveniently for e -» or 00. Let us assume for e — > 

^i(e, x) -> (137) 
- ^[(/o + /ilne + / 2 ln 2 e ) + e(^ + ^ln e + ^ln 2 e) + ...] , 

where p can be integer or half integer and f g , g q , . . . depend on the particular 
process. For the various jet shape variables in e + e~ annihilation (thrust, jet 
mass, C parameter) e.g., .Fi(e, x) — > -^fy^/e with fy = 4, 2, 6ir respec- 
tively. 

7 Note that in the original papers the quantity <5sgd(w 2 ) was denoted as a e ff(w 2 ) . 
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As a consequence we have for Q — ► oo 

[ «) - /„£ + h - £ b?) + . . . (138) 

or 

y(Q 2 ,x) = ^ PT (Q 2 ,:r) + ^[CiA 2p + C 2 A' 2p + C 3 A'> p ] (139) 
where we have introduced the moments 

A2p " 2vr 7 m 2 m sgdI J , 

4* = T- r Wfi§g D (m 2 ) , (140) 

Z7T Jq 771 

«/ o 

The coefficients Ck are dimensionless, are in practice at most linear in In Q 2 
and are calculable but process dependent, A 2p , A' 2p , . . . are theoretical un- 
known but should be universal. They could be determined in principle by 
studying the dependence on the hard scale Q of appropriate observables, 
V(Q 2 ,x). 

Terms in 1/Q should emerge for what we have seen in various e + e~ jet shape 
variables; correction of the type 1/Q 2 in the DIS structure function (see 
[106J). The results can be equivalently expressed in terms of the moments 
of <ysQ D (Q 2 ) rather then of a^ G (m 2 ). The experimental situation is not 
completely clear but the data seem to be consistent with 

Al = T- r IW® c4 P G (Q 2 ) = -A x « 0.2 - 0.25 GeV (141) 
2ir J Q Q 2 7i 

and 

A2 = T- r ^Q^sEciQ 2 ) = -A' 2 - 0.2 GeV 2 . (142) 
^ Jo Q 

Note that in principle the quantity Vpx(<3 2 , x) is well defined, since «sdg has 
no unphysical singularities, due to its definition (|S9*jl. To be consistent up 
two-loop one should use equations for a^ G of the type considered in sec. 4.1 
or 4.2 with the appropriate value of Asgd- However, often an infrared cutoff 
Qi has been introduced in the application and A\ related to the quantity 
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( a ) considered in sec. 3.1. 

Relations involving OPE and the non trivial IR structure of the theory have 
been considered in the 2D Gross- Neveu model in |lU7j . 
To conclude this section we should mention that the same physical problem 
we have briefly discussed from the point of view of the running coupling is 
the object of another very alive line of research that operates in the frame- 
work of the Borel summability and the renormalon singularities. Since there 
are obviously connections between the two perspectives, we must invite the 
interested reader to consult some of the existing excellent general reviews, 
e.g. |l08j . and among the most recent works we remind ref.[109 . 

4.6 Analytic Pertubation Theory 

As explained at lenght, analytization displays a suitable method to get rid 
of unphysical singularities which affect the standard expansion parameter 
in t- and s-channel. Thus the issue of how perturbation theory should be 
accordingly modified naturally rises, and it has been investigated from phe- 
nomenological and theoretical point of view (e.g. |110j ). We briefly recall 
here a consistent way to incorporate the ghost-free model for the IR finite 
couplings eqs. (jl()5j) and ()125j) within perturbation theory, known as Analytic 
Perturbation Theory (APT), and developed in |lllj - [TT4*j . The main require- 
ment here is the subtraction of unphysical singularities in the RG improved 
series for physical observables as a whole, by computing their discontinu- 
ity across the time-like cut, as it has been done for the space-like coupling 
itself eqs. (jl()4j) and (jl()5j) . Specifically |113j . given a space-like observable 
perturbatively known 




(143) 



k>l 



one can define the k-ih spectral density 




(144) 



Then eq. ()143J) is to be substituted by the ghost-free expansion 




(145) 



k>l 



where, leaving the loop level understood, 




(146) 
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Then a standard power expansion is converted into a non power one. Clearly, 
there is now no unique expansion parameter, but an entire set of ghost-free 
expansion functions (j!46|) at any loop level, each defined by the analytization 
of subsequent powers of the perturbative coupling. Obviously from (j!46|) with 
k = 1 the analytic coupling (jl05|) is recovered at each loop level. This recipe, 
working for space-like observables, is manifestly quite analogous to the RKP 
non-power expansion (J56j) of a time-like observable; this in turn, within this 
framework, can be rexpressed via the k-th spectral density (j!44|) [113 

R APT (s) = l + J2 d Ms) (147) 

fe>i 

and 

0*00 = -/ — Pk(a) (148) 

which yields for k = 1 the analytic time-like coupling itself ()125|) . The key 
point here |113j is that, due to the forced analyticity of the coupling and its 
analytized powers f)146|) . the two sets ()146|) and ()148|) are put into one-to-one 
relation by the linear integral transformations (J53|) and (|54j) . in this context 
usually renamed R and D respectively, that is 

A k (Q 2 ) = D[H*(s)] , 51*00 = R[A(Q 2 )] • (149) 

This yields a closed theoretical scheme for representing observables of any 
real argument, both space-like and time- like (for a quite recent review see 
|114j to which we also refer for technical details). The main features of 
these functional sets are illustrated in fig^ taken from |1 14j . As yet noted 
the first function in both sets coincides with the relative analytic coupling, 
respectively ()105|) and ()125|) . while Ak>2 and 2lfc>2 start with an IR zero 
and oscillate in the IR domain around k — 1 zeros; furthermore they all obey 
the UV asymptotic 1/ ln fc z resembling the corresponding powers ce^(z) |114j . 
The differential recursion relations [89j 

idA^m_ ^ (n) i gW(g) ^ (n) 

j=i i=i 

(150) 

that hold in all orders jHI], allow to relate different analytized powers within 
each set, tough explicit expressions are reliable in a simple form only in the 



L 




Figure 4: (a) Space-like and time-like global analytic couplings in a few GeV domain 
with rif — 3 and A = 350 MeV; (b) "Distorted mirror symmetry" for global expansion 
functions, corresponding to exact two-loop solutions. 

one-loop case fe.g.[8*9"]). Meanwhile, Ak and 21^ for k = 1,2,3 have been 
tabulated in jHO] up to three-loop, starting from the exact two-loop solution 
(J23|, and the Pade' approximant technique in the three- loop case; further 
improved approximations as eq. (JH!Z|) have been exploited at low scales in jHI] 
(see also |119j ). These numerical values are reproduced quite well in the 
range (2, 100) GeV by "one-loop" inspired approximate formulas, derived for 
practical aims in |89j, and depending on an effective scaling constant A e ^ 
as a fitting parameter (see Tab.l in ref. |89|). Computation of the functional 
sets has been recently extended to include analytic images of any real power 
of the coupling (Fractional Analytic Perturbation Theory) on the grounds of 
properties of the transcendental Lerch function |120j . 

The APT algorithm can also include thresholds effects |113j in a real de- 
scription, by modifying the k-th spectral density ()144|) discontinuously at 
the heavy quark thresholds irtf 

Pk{<?) = Pk(cr, 3) + 6 (a - rrij) [p k (a,n f ) - pk(cr,n f - 1)] , (151) 

n/>4 

descending from the trivial matching condition (cf. sec. 4.2). The global 
functions resulting from densities (|151|) . Ak arid 21^ , can be obtained from 
the local ones with rif fixed, by adding specific shift constants Ck{nj) , not 
negligible in the rif = 3,4 region (see |113j for numerical evaluations; for 
instance, in both t- and s-channel it has been estimated ci(3) ~ 0.02 ). 
The main tests of APT being obviously at low and intermediate scales, 
a number of applications to specific processes, both in the space-like and 
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time-like (low and high energy) domain, have been quite recently performed. 
For instance, Bjorken and Gross-Llywellin-Smith sum rules r-lepton 
|112| (see also |115j ) and Ypsilon decay [55] . e + e~-annihilation into hadrons 
[53] IH2l I114j and hadronic form factors [1161 1117j . As a result (see tab.l in 
ref. [118j ). the main advantages of the APT approach are better convergence 
properties of the ghost-free non-power expansion than the usual power one, 
since the three-loop term is always strongly suppressed (even less than data 
errors). This feature yields a reduced scheme and loop-dependence. Finally, 
transition from Euclidean (space-like momentum) to distance picture has 
been also developed in |121j to which we refer for details, involving a suit- 
able modified sine-Fourier transformation, consistently with the APT logic. 

5 Lattice Theory 

In all development above one starts from manipulation on the perturbative 
expansion and try to remedy to its lack of convergence, but one can never 
attain to really non perturbative effects related to the singularity of the the- 
ory in a s = 0. 

Up to now, the only general technique to handle (even if with many limita- 
tions) non perturbative problems is Lattice Theory. 

We want here to mention essential ideas and summarize the present status 
of art for what concerns the running coupling. 

We refer to the standard formulation as due to Wilson |122j . a classical text- 
book is e.g. |123j . 

On the Euclidean lattice the spacing a plays the role of the UV cutoff of the 
continuum formulations. In momenta space the fields are then defined on 
the Brillouen zone but, for practical purposes, the lattice is finite, thus an 
IR cutoff appears too and the momenta become finite. The quantum theory 
is obtained via path integral quantization and the result is similar to a sta- 
tistical mechanics formulation in D = 4, at an inverse temperature (3 which 
turns out to be 2N/g$, for the SU(N) gauge theory and where go is the bare 
coupling constant. The lattice is a gauge invariant regulator but for finite 
lattice the functional integral is perfectly convergent, and no gauge-fixing 
procedure is in principle necessary to evaluate gauge invariant quantities. 
Very used are improved actions, namely actions with extra terms of physi- 
cal dimensions higher than four, with coefficients suitable chosen in order to 
have a better approach to the continuum limit (see |124j for a review). A 
complication of the lattice formulation is the well known doubling problem 
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of the fermions, in order to have in the continuum limit the correct number 
of fermions one has to use an action which breaks the chiral symmetry for 
finite a. The chiral symmetry can be restored by considering a critical value 
for the hopping parameter for any lattice spacing. Recipes for the fermionic 
action which differ from the Wilsonian one are widely used, each one with its 
own problems. To get answers from the lattice usually one does not compute 
the functional integral but perform Monte-Carlo simulations. The compu- 
tation of the running coupling constant on the lattice starts from its non 
perturbative definition in the low energy region; there are a lot of definitions 
which have been actually used, and in this region they can have a very dif- 
ferent behavior with the energy scale, but extrapolated in some way at high 
energy they must behave similarly, and it has to be possible to correlate one 
to another using perturbation theory. That this happens is a consistency 
check of the theory on all the energy scales. In particular at high energy it is 
possible to commute to the MS scheme in dimensional regularization, defined 
only in perturbation theory. An exhaustive review on the various possible 
approaches is in [125] . 

From a physical point of view the framework of the computation is fixed by 
setting the scale. In the pure gauge theory one needs, from the experiments, 
the knowledge of a physical quantity of mass dimension Mh, and on the 
lattice, starting from a reasonable go, or equivalently f3, one has to be able 
to compute the number aM H , and this determines the spacing a in physical 
units. The window of energies on the lattice, in which M# has to fall is then 
[1/L, I /a]. L/a is the number of sites for side of the lattice and so is fixed 
by the computer power, typical values are 16 — 64. 

The asymptotic freedom can be formulated also in terms of bare coupling 
constant which vanishes for a — > 0, for instance we have the two loop for- 
mula 



(the ratios A^/A, where A refers to a suitable continuum scheme, e.g. MS, 
have been evaluated in lattice perturbation theory long time ago |126j . |127| ). 
Being interested in the continuum limit one would like to have a and go as 
small as possible but one needs to take into account the limitations of the 
window energy. The dependence on the chosen observable should be not 
important as far as we are "close" to the continuum limit (scaling region), 
used quantities are the string tension, the mass splitting in heavy quarko- 
nium or the hadronic radius r |128j . In the nineties there have been several 
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determinations of the running coupling constant in the pure gauge theory or 
in quenched QCD ([T2g]-fnffi|). 

Among the papers, |135| and 136] are closer in the formalism to the other 
parts of this review, indeed the definition of the renormalized coupling con- 
stant uses a condition on the trigluon vertex, therefore very similar to that 
employed in usual perturbative procedures. In principle such a definition of 
the coupling can however be used in the full infrared region, the only lim- 
itations being due to the lattice size effects. The disadvantages of such an 
approach are that A p on the lattice is a rather unnatural quantity, defined in 
conventional way through the fundamental link gauge variables; moreover, to 
compute the Green functions, the gauge fixing is needed. For the implemen- 
tation of the gauge fixing on the lattice a good review is |137| . It is worth to 
mention that even the lattice formulation of the Landau condition suffers for 
the problem of the Gribov ambiguity. Actually one could pick out a unique 
element in each gauge orbit by finding the absolute minimum of a suitable 
functional with respect to gauge trasformations but, from a numerical point 
of view, to distinguish the absolute minimum for instance from the local 
ones is a very difficult task. Beyond |137j a clear analysis of the numerical 
situation is in |138j . relations with other approaches are in the review [139J, 
devoted to the infrared behavior of the QCD Green functions. On the lattice, 
by Monte Carlo average, after performing the Fourier transform and in the 
Landau gauge, one computes the unrenormalized Green functions 

4 n L.M>l- >Pn) = <4*(Pi)- >4*(Pn)>- ( 153 ) 

In |135j and |136j the scale is set by a previous determination of the string 
tension |129j . then /3 = 6 and correspondingly a" 1 = 1.9 ± 0.1 GeV in a 
hypercubic lattice are taken in |135| . more lattices are used in [136J. 
Assuming that finite-volume effects and discretization errors are under con- 
trol, for suitable small momenta one can adopt the formalism of the contin- 
uum QCD and then to write for the gluonic propagator 

GuM=Gu( P 2 )^-^) (154) 

and for the three point function (after factorization of the color tensor f a b c ) 
Gu\ up {Pl,P2,P3) = ^ ( u a( 3 7 {Pl,P2,P3)Gu a ^{pi)Gu f3u {p 2 )G u lp {p 3 ) . (155) 

Where Pi + P2 + P3 = 0. On the lattice of course one measures the Green 
functions, so differently from the usual perturbative approach the renormal- 
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ization conditions are expressed on the Green functions and not on the 1-P-I 
functions. Thus 

G R (p)\ p 2 =fl 2 = Z g 1 (a^)G u (pa)\ p 2 =l ,2 = — (156) 

is the condition which fixes the wave-function renormalization constant in 
non perturbative way. For the vertex the most obvious approach would be 
to impose a condition on the most symmetric point (MOM scheme), but 
from a numerical point of view it turns out that an asymmetrical condition 
in which one of the momenta vanish works better (MOM in continuum QCD 
|126j ). Defining the vertex renormalization constant according to 

Gu (pay Gu (o) = 6 ^ {pa)9 ^ ' (157) 

the vertex renormalization condition reads 

g(n) = Z 3 J 2 (an)Zy 1 (an)g . (158) 

In a massless theory the Green functions at exceptional momenta are po- 
tentially divergent, but in this case all is finite |14()j . [T4"T] . Note that the 
tensorial structure of the previous formulae is not obvious because on the 
lattice 0(4) symmetry no longer holds; a test of consistency of their numeri- 
cal computation performed by the authors is to check the previous structures 
for some values of pa . As a matter of fact the parametrization of the for- 
mulae (J154J) and (|155|) cannot be the most suitable in order to extract the 
coupling, lattice perturbation theory would suggest using p il = \ sin(^) in- 
stead of p M , indeed the authors from |136j claim that only in this case they 
are able to find evidence of the running of the coupling. In order to prove 
the running of the coupling one fits the data of the simulation with the curve 
which express A through g(n), see FigJSJ using the two loop formula (but in 
|13bJ they control the three loop effect) 

l -A 

A m5m = /^exp(- 2^^) )^°^)) ^ ( 159 ) 

in both the papers one finds an evidence of plateau close to 2 GeV, see 
FiglSfb), beyond A(g) decreases due to cutoff effect. The effects of finite 
size are controlled by performing the computation using more lattices with 
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Figure 5: (a)The QCD coupling constant a M Q M from |136| . \i is in GeV. The full line is 
the three-loop running. (b)A from |136| . 

the same lattice spacing. 

Having computed A M g M one can get Aj^g. To this end one has to do a 
one-loop [213 calculation which expresses perturbatively the relation between 
dMsif 1 ) an d #mom(/-0- Indeed a discrete renormalization group transforma- 
tion which connects the Green functions computed on the lattice and those 
computed in MS scheme must exist. Being this scheme on the lattice specified 
by conditions on the Green functions, the transformation is easily determined 
by the one loop MS continuum calculation of the trigluon vertex and of the 
propagator at the MOM point and in Landau gauge. More difficult compu- 
tations in lattice perturbative theory are not needed in this approach. 
In |135j it is obtained A^jg = 0.34 ± 0.05 GeV; the authors of [136 , with a 
writing which distinguishes the possible errors coming from the lattice spac- 
ing, instead get A^ = (0.303 ± 0.05 GeV) 1 ^ y . 

The same renormalization scheme has been applied to measure the coupling 
using the gluon-quark vertex in quenched simulations in |142j . 
In order to avoid the major source of systematic errors, in spite of the com- 
puter cost, simulations with dynamical quarks have been addressed. The task 
has been undertaken using NRQCD |143j in the past; also the authors of 136 
published a computation with nf = 2 dynamical quarks |144j in the MOM 
scheme previously described. Interestingly the matching with the perturba- 
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tive behavior of the coupling, and then the extraction of A, was made taking 
into account an effect of OPE condensate arising from the gluonic propagator; 
therefore the data computed on the lattice for the coupling where compared 
not with a S)Per t(^) but with a S)Per t(j)(l + -%), where c has to be adjusted in 
a combined fit. The form of this correction, similar to those introduced in 
eq. (|128|) . is based on previous lattice computations and already utilized in 
order to improve the determination of a s in quenched simulation |145j . 
They also perform the extrapolation of the coupling to nf = 3 and then the 
renormalization group evolution to the value of Mz of the energy through 
the thresholds, FigEfa) quotes their result. 

A recent review |146j analyzes the results from lattice calculations up to 
2004 with emphasis to the theoretical issues in dynamical quark simulations. 
In particular a comparison among determinations of the various coupling 
constant is presented, besides the already cited papers, FigJHfa) includes 
some newer |147j . |148j and some older result ,150j. The most recent re- 
sult in FigEta), obtained by the HPQCD collaboration (2004), is a s (M z ) = 
0.1175(15) in good agreement with the PDG 2004 average 0.1182(20) [IE]. 
It has be noted that the results obtained using Wilson fermions |148j . |149j 
are significantly lower than those obtained using staggered fermions |147j . 
An original approach to the lattice computations has been developed by 
the ALPHA collaboration which recently performed a determination of the 
strong coupling constant with two dynamical flavors |151j . Exposition of the 
method is in |152j , a detailed discussion of the mass renormalization in |153j . 
One of the computational problems in the usual approaches is that the win- 
dow energy must be wide to contain both the matching to the hadronic scale 
and the asymptotic behavior of the coupling. In the approach of the ALPHA 
collaboration one can determine a recurrence relation between quantities re- 
ferring to double energy each step, and thus one constructs a ladder from the 
hadronic to the perturbative regime. 

The technical tool is an intermediate renormalization scheme in which the 
renormalized coupling constant is defined by derivative with respect to pa- 
rameters characterizing the boundary conditions of the Schroedinger func- 
tional, namely the amplitude for a transition from a given field configuration 
at a given time, t — 0, to another one at t — T. Moreover the bound- 
ary conditions yield a natural infrared cut-off so it is possible to perform 
simulations at vanishing quark mass |153j (using Wilson fermions the mass 
renormalization is additive, and imposed using the PCAC relation). The 
renormalization of the coupling is defined on the massless theory thus a non 
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perturbative mass independent scheme is obtained. In other computations 
the light quark masses are actually rather large, and chiral extrapolations 
are needed in order to compute the coupling to be compared to the MS one. 
Taking T = L, being a and L the only dimensional quantities, the following 
relation between bare and renormalized couplings holds g L = g(a/L,g ) . 
Here 1/L plays the role of the finite scale of energy. The same procedure 
can be applied by doubling the size, L — > 2L, and after elimination of go 
one finds a relation between the two renormalized coupling constants. At 
the end, by extrapolating the limit a/L — > 0, they define the "step scaling 
function" which connects the two couplings: 

(?{u L /2) = u L . (160) 

In f!160|) u = g 2 and the two couplings refer to lattices of size L/2 and L at 
fixed Physics, namely bare parameters. The previous equation is actually a 
discrete renormalization group transformation, indeed a relation between the 
j3 and the a functions is determined. Of course the step scaling function is 
constructed by a numerical procedure; consistency with perturbation theory 
requires a{u) ~ u as u — > 0, for larger u, a{u) is computed by a suitable 
interpolation starting from the set of points determined on the lattice (a 
discussion of the systematic error arising from this interpolation, which it is 
claimed to be negligible in a suitable range of coupling, is in |154j ). 
Starting with a given coupling and a given size (u max , L max ) using the step 
scaling function one computes Ui = g 2 (2~ l L max ) and for sufficient high en- 
ergy, with the identification fi = 2 l /L max , one matches the pertubative regime 
evaluating A (the ratio between this A and Ajyjg is known). 
The starting point of the iteration is a priori unknown, it is determined 
by matching some hadronic quantity, as a practical approach it was com- 
puted r in the chiral limit; assuming r = 0.5 fm it turns out A^r 2 = 
245(16) (16) MeV. For the low energy behavior of this coupling see [157 . 

The most recent determination of the running coupling constant has been 
obtained by a HPQCD — UKQCD collaboration |155j . Using a new dis- 
cretization of the light quark action the authors perform the first computa- 
tion with three dynamical quarks. After having fixed the lattice parameters 
(m u = rrid, m s , m c , rrib and a) matching five experimental quantities the 
determination of the coupling is obtained by an optimization of lattice per- 
turbation theory, in which the coupling must fit 28 short-distance quantities 
calculated on the lattice. The final estimate is a s (Mz) = 0.1170(12), which 
slightly corrects the already reported result of |147j . 
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Figure 6: (a) Strong coupling constant a s (Mz). Lattice data are from Aoki et 
al.[ISI|, Davies et al.QI3], SESAM|T23|. Boucaud et al.[ffi|, QCDSF-UKQCD|T1%| and 
its updated]], HPQCD-UKQCD-MILC-Fermilab^^ and its update EH] . (b)Values 
for av versus d/a from each short-distance quantity from |155j . with (top) and without 
(bottom) light-quark vacuum polarization. The dashed lines show the predictions of the 
perturbative running which best fit the distributions. 



The figure E(b), from the cited paper, shows the remarkable effect of the vac- 
uum polarization of the three light quarks on the coupling ay introduced in 
eq . ([67)1 (the momentum is represented in units of the inverse lattice spacing). 
Among the many lattice results, we would like now to pay a little attention to 
those concerning the behavior of the coupling in the very low energy region. 
The MOM definition of the coupling is likely the best suited to this aim, 
indeed the already cited authors using this approach addressed this task in 
|15b'j . In order to reach low momenta also lattices with a rather large lattice 
spacing, and then low (3, have been used. The results for the pure gauge the- 
ory are shown in Fig|7| The number of data in figJTfb) is sufficient to check 
the theoretical model of the vacuum as instanton liquid. According to this 
model the coupling must have a p A dependence for small energy, the data in 
FigEfa) are indeed fitted by this power low, the proportionality coefficient 
yields the instanton density n: 

<*mom(p) = T5~ n- V- (161) 

187T 

Such a behavior is not a universal feature, from a glance at the figures 8 
and 9 in the second paper of |142j . in which the coupling is defined in MOM 
scheme but using the quark-gluon vertex, one realizes that there the coupling 
is still vanishing for small momenta but, in spite of the rather little number 
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Figure 7: (b) Symmetric MOM coupling constant for different lattice settings and fits 
to perturbative expression plus power corrections in the high momenta region and to the 
expression of ea. l|161(l for small momenta, (a) Region of small momenta is zoomed. 

of data, it does not seem with the p 4 low. On the contrary the coupling 
cey(Q) introduced in eq.flHZJ) should instead diverge for Q — > 0. 
In the last years infrared features of QCD has been addressed using the 
Dyson-Schwinger equations. This approach, which starts from suitable trun- 
cations of the DSE, seems in particular effective in the study of the low energy 
behavior of the propagators (see |139j for a review). In this way infrared ex- 
ponents for the gluon-gluon and ghost-ghost two point functions have been 
determined. In the Landau gauge using the ghost and the gluon propaga- 
tors one can give a definition of the invariant coupling (which then resembles 
the familiar definition of QED) and it turns out that it has a finite infrared 
limit |158j . It is found evidence for this fixed point at least in the context 
of pure SU(2) lattice gauge theory |159j . however this conclusion does not 
seem supported in other recent lattice SU(3) simulations 160j- [T6*T] . The two 
results are compared and discussed in |162j . very recent results also in [163J. 
As yet noted, such a lack of universality of the running coupling constant 
in the low energy regime is somewhat expected, however on such different 
behaviors, which are a little disturbing, some comment is in order. In |164j 
a comparison of various approaches is made with reference to singular RG 
trasformations. The RG transformations are associated to a change of vari- 
able g — ► g' = f(g), which in general depends on the dimensional parameters 
too; in the perturbative framework f(g) is usually supposed to be an analytic 
function in the origin g = 0. Now considering the MOM schemes, it is well 
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known that in perturbation theory the definitions of the coupling are vertex 
dependent |2Hj but, at least in the massless case, this dependence turns out 
to be weak, namely the function / is close to the identity It is conceivable 
that in the infrared region, all more reason taking into account the different 
dependence of the coupling on the quark light masses, the function / can 
have non trivial consequences. 

Some proposals that we have examined which tries to extend the perturbative 
series in the low energy region can be seen as a not trivial RG transforma- 
tion, at the beginning defined in full perturbative region, but such that the 
perturbative series in the new variable does not suffer from some pathology 
of the first one. Thus formula ([11)7)1 . which avoids the Landau pole, can be 
viewed, as already noticed, as the result of the non-analytic transformation 
a s -> a' s = a s + ^(1 - exp(^))- 1 . 

Every devisable transformation should share with the previous one the two 
nice properties to enjoy the asymptotic freedom, a s ~ + — > a' s ~ + , and to 
be free of ghost singularity, a s = oo — > a' s = finite. A priori there are not 
other compelling requirements, and one can imagine a situation in which, 
with respect to a given problem, a definition, and then a consequent pertur- 
bative expansion, has better properties than another. The considerations in 
|164j show that with other transformations one can induce various infrared 
behaviors in the new coupling. The conclusion is that "there is not direct 
physical sense in attempts to establish some "correct IR behavior" of the 
perturbative QCD invariant coupling". 

6 Comparison with the data and conclusive 
remarks 

For completeness in this concluding section we want to comment briefly on 
the experimental situation as continuously kept updated by jTH] and other 
specialized reviews (see e.g. 0). Tab.l in the contribution of S.Bethke in 
this issue summarizes the determination of a s (Q 2 ) at various energies be- 
tween 1.58 and 206 GeV extracted by various types of experiments: deep 
inelastic scattering (DIS), e + e~ annihilation (hadron cross section, structure 
function, jets and jet shapes), hadron-hadron collision, Z, T and r decay, 
heavy quark bound states. 

In fig. |Bf a) from jH] such results are compared with the 4-loop perturbative 
expression of cty^(Q 2 ) as given by (|H1|) together with threshold matching. 
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Figure 8: (a) From summary of measurements of a s (Q 2 ) . (b) From compari- 
son between the 1-loop analytic coupling i|107|) with AS nf ~ 3 ^ = 206 MeV, and the values of 
a s fitting experimental data on quarkonium spectrum within Salpeter formalism (circles, 
pentagrams and squares refer to light-light states, diamonds and crosses to heavy-heavy, 
plus signes and asterisks to light-heavy states). 



As it can be seen the agreement is very good inside the errors and this con- 
stitutes a very significant test for QCD. Note, however, that assuming all 
expressions equally normalized at Mz = 91.2 GeV the difference between the 
3-loop and the 4-loop expressions is of order 1/10000 in the interval from 10 
to 200 GeV, it becomes of the order 1/1000 between 2- loop and 4-loop and of 
few per cent between 1-loop and 4-loop. We have in this way the proof that 
at three loops the theory is practically at convergence in comparison with 
the precision of the experimental data and that two loops (naturally after 
the appropriate rescaling of A) is already a very good approximation in the 
considered range 8 . The choice of a particular RS or another is essentially 
immaterial. However, under a scale of few A^ n f =3 ^ (let us say 2 GeV), the 
MS becomes useless, due to the Landau singularities, and we have to refer to 
any alternative scheme free of unphysical singularities and to follow it very 
consistently. A comparison with the data becomes RS dependent and would 
have no meaning out of the well defined framework. 



Obviously the difference explodes near the singularity. 
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Note that the analytic scheme discussed in sec. 4 seems particularly conve- 
nient at this aim. It shares the simplicity and the universality of the MS 
scheme at a large extent. It provides a coupling which is regular in the entire 
interval < Q 2 < oo and has a finite limit for Q — ► independent of the 
number of loops used. As a consequence of this fact its expression converges 
much faster as the number of loops increases. The difference between the 
2-loop and the 3-loop expression is again of order 1/1000 but this time more 
or less in the entire range from to +oo. Even the difference between the 
1-loop and the 3-loop is of order of 1/100 and may be a fraction of this in 
limited ranges. 

For construction a an (Q 2 ) coincides asymptotically with oiy^{Q 2 ) for large Q 
and ain (Q 2 ) fits the data of Tab.l from nearly as well as a^(Q 2 ). Note 
however that, as seen in sec. 4.2, we have more precisely 



The coefficients of seem to decrease as the number p of the loops increases. 
However for the value of p we actually use the power term in ()162|) can be 
non negligible at some intermediate scale and may be disturbing in particu- 
lar with reference to the problematic discussed in sec 4.5. As we have seen 
they may be eliminated for every given p by choosing an appropriate non 
perturbative term a ai f (Q 2 ) to be added to the "perturbative" o^{Q 2 ). As 
we have seen phenomenologically modified expressions for a an (Q 2 ) have also 
been proposed to eliminate such term. 

Information on the QCD interaction under few GeV can be obtained from 
power corrections to various observable at some intermediate energy, as we 
have seen. Other information can come from relativistic calculation of the 
spectrum of light-light, light-heavy and from the highly excited heavy-heavy 
quarkonium states. 

In |165j . e.g., a reasonable reproduction of the entire calculable spectrum and 
in particular of tt — p, the K — K* and r] c — J /ip was obtained. Such analysis 
was performed in the framework of a Salpeter formalism constructed using as 
input only a generalization of the ansatz (|64j) and a 1-loop analytic coupling 
with n f = 3 and A( n /= 3 ) = 180 MeV. Note that such a£2(Q 2 ) differs from 
ot£2(Q 2 ) with A(™/= 3 ) = 375 MeV by less than 0.5 %. 

Actually the developments reported in sec 4.5 were made in the framework 
of the SGD scheme and this scheme may be the most natural even for bound 
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states problems. In fact the inclusion of certain second order corrections 
would simply amount to a rescaling of A. However the SGD and the analytic 
scheme are strictly related and have even a common root in a sense |105j. 
They are both based on dispersive techniques and can be made to coincide 
up to two loops. 

In any way according to a more recent preliminary analysis, made in the 
context of the mentioned relativistic determination of the quarkonium spec- 
trum |166| . it seems that the curve ai}(Q 2 ) (with A^ =3) = 206MeV) is in 
very good agreement with the data up to 200 MeV corresponding e.g. to the 
cc(lD) state (see figlEtb)). On the contrary for higher excited states, corre- 
sponding to 100 < Q < 200MeV, the values of a s that reproduce the data 
seem to be somewhat lower than ain (Q 2 ) and this could confirm the need of 
a non perturbative contribution to be added to a an (Q 2 ) as given in sec. 4.5 
[166J. However, the experimental situation is so uncomplete and confused 
in this range and the applicability of the theoretical method so questionable 
that no real significant conclusion can be drawn. 

In conclusion, the use in QCD of a RS, in which the running coupling con- 
stant is free from Landau singularities, seems to offer a framework in which 
all existing phenomenology can be discussed and non perturbative correc- 
tions usefully parametrized. It is even a framework in which a comparison 
with Lattice Theory results is easier and more natural. 
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